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Preface 

This  investigation  was  sponsored  by  the  Flight  Control  division 
of  the  Air  Force  Flight  Dynamics  Laboratory,  Wright-Patterson  AFB. 

My  thanks  are  expressed  to  Mr.  M.  A.  Ostgaard  of  that  Division  for 
his  helpful  discussions  on  the  problems  of  flight  control,  and  for 
his  assistance  in  providing  me  with  an  introduction  to  current 
literature  on  the  subject. 

In  addition,  I  want  to  thank  Mr.  F.  J.  Sansom  of  the  Director¬ 
ate  of  Computation  for  his  advice  on  the  use  of  the  MIMIC  programm¬ 
ing  system  in  the  IBM  7094  computer,  and  also  my  thesis  advisor, 
Professor  C.  H.  Houpis  of  the  Air  Force  Institute  of  Technology, 
for  his  support  and  encouragement. 

It  has  been  said1  that, 

A  decade  of  promise  and  a  dearth  of  actual  applications 
mark  the  turbulent  past  of  self-adaptive  flight  control. 

Presented  in  this  thesis  is  an  outline  of  the  problem  involved 
in  designing  a  control  system  for  a  high  performance  aircraft  or 
aerospace  vehicle,  and  an  indication  is  made  of  one  possible 
avenue  of  ajjproach  to  its  solution  which  may,  perhaps,  be  worth 
further  development. 

For  the  student  of  control  theory  or  electrical  engineering 
who  iB  unfamiliar  witli  the  particular  field  of  aircraft  control 
systems,  the  article  referenced  below  is  recommended  as  an 
introduction  to  the  subject. 


1Andeen,  H.  E.  "Self-Adaptive  Autopilots."  Space/Aeronautics  1 
pp.  46-52,  April  1965. 
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Abstract 


In  a  high  performance  airoraft  ohanges  in  Vaoh  number ,  angle 
of  attack  and  altitude  can  oause  a  large  variation  in  the  short- 
period  transfer  function.  To  provide  the  pilot  with  a  constant 
pitch  rate  control  characteristic  an  airborne  computer,  with 
inputs  of  elevator  deflection  angle  and  pitoh  rate,  is  UBed  to 
identify  and  track  changes  in  the  elevator  effectiveness. 

Empirical  equations  are  defined  to  approximate  the  airoraft 
time  constant,  damping  factor  and  natural  frequency  as  functions 
of  elevator  effectiveness  in  three  difference  equations,  which 
are  iterated  to  model  the  airoraft*  Parameters  in  the  difference 
equations  are  perturbed  until  the  equation  which  uses  a  value  of 
elevator  effectiveness  intermediate  between  the  values  in  the 
other  two  equations  also  has  the  smallest  mean  square  error  from 
the  aotual  aircraft  response.  The  value  of  elevator  effectiveness 
in  this  intermediate  equation  is  then  presumed  to  be  the  same  as 
that  of  the  aircraxt,  and  is  used  to  set  the  loop  gain  to  a  pre¬ 
determined  suitable  value. 

Simulation  with  an  aiZ'Oraft  whose  elevator  effectiveness 
varied  over  a  range  of  240 il  showed  that  the  desired  loop  gain 
was  maintained  within  a  faotor  of  two  for  both  pilot  command 
inputs  and  for  random  wind  gust  disturbances  of  root-mean-square 
amplitude  20  ft. /see.  Gain  adjustment  rates  of  up  to  6  dB/seo. 
were  achieved. 
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A  SELF-ADAPTIVE  AIKCRAFT  PITCH  RATE  CONTROL  SYSTEM  EMPLOYING 
DIFFERENCE  EQUATIONS  FOR  PARAMETER  IDENTIFICATION 

I.  Introduction 

Feedback  Design  Philosophy 

For  a  plant  which  has  well  defined  and  relatively  constant 
parameters  the  techniques  for  designing  a  closed  loop  control 
system  which  has  a  fast  response  and  adequate  stability  have  been 
extensively  formulated.  These  techniques  ensure  that  small  amounts 
of  drift  in  the  values  of  the  plant  parameters,  due  to  ageing  or 
environmental  changes,  have  an  insignificant  effect  on  the  overall 
performance. 

However,  in  the  case  where  the  gain  of  the  plant  (as  distinct 
from  its  dynamic  terms)  has  a  wide  range  of  variation,  the  design 
problem  becomes  aoutei  if  the  loop  gain  is  set  to  a  value  which 
provides  satisfactory  operation  when  the  plant  gain  is  high,  the 
control  response  will  be  slow  when  the  plant  gain  decreases;  but 
if  the  loop  gain  is  set  to  a  value  suitable  for  low  plant  gains,  the 
system  may  become  unstable  when  the  plant  gain  is  high.  One  way  of 
reconciling  these  conflicting  requirements  is  to  use  lead  compen¬ 
sation  to  increase  the  frequency  at  which  the  loop  phase  shift 
readies  -l8o  deg.  When  suoh  compensation  is  made,  it  is  theoretic¬ 
ally  possible  to  aohieve  the  desired  loop  gain  for  low  plant  gains, 
with  sufficient  gain  margin  to  accommodate  the  maximum  plant  gain. 

In  principle,  for  a  truly  linear  control  system,  any  degree  of 
plant  gain  variation  may  be  tolerated  if  enough  lead  compensation  is 
used.  Practical  considerations,  however,  impose  restrictions  on  the 
effective  amount  of  compensation  which  can  be  applied.  As  the  plant 
gain  increases  from  its  minimum  value,  so  alec  does  the  control 
system  bandwidth,  so  that  for  large  plant  gains  the  bandwidth  is 
greater  than  necessary.  The  result  is  that  the  oontrol  system  has 
undue  sensitivity  to  noise  arising  from  spurious  inputs  to  the  loop, 
and  from  the  feedback  sensors.  Additionally,  high  order  plant 
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dynamios  cause  phase  shifts  at  high  frequencies  which  tend  to 
destabilize  the  loop.  Thus,  in  the  event  of  large  variations  in 
the  magnitudes  of  the  plant  parameters,  an  alternative  approach  to 
control  system  design  is  indicated. 

Adaptive  Control  Systems 

A  high  performance  aircraft  is  an  example  of  a  plant  for  which 
the  design  of  a  control  system  entirely  by  conventional  means  is 
difficult,  if  not  sometimes  impossible,  without  degradation  of  the 
handling  qualities  in  some  flight  regimes.  This  situation  arises 
largely  because  of  the  variation  in  control  surface  effectiveness 
as  a  result  of  changes  in  Mach  number,  air  density  and  angle  of 
attack,  etc.  Design  of  a  control  system  in  which  gain  is  scheduled 
as  a  function  of  the  output  of  an  air  data  computer  is  one  solution. 
A  system  of  this  type  is  called  adaptive,  since  it  is  continuously 
altered  automatically  to  suit  the  surrounding  environment.  The 
major  limitation  is  that  the  relationship  between  aircraft  response 
characteristics  and  particular  Plight  Conditions  (P.C.)  is  indirect, 
so  that  a  completely  satisfactory  gain  scheduling  program  can  only 
be  arrived  at  after  extensive  flight  testing. 

Self-adaptive  systems,  in  whioh  gain  (and  possibly  other 
control  parameter)  adjustment  is  made  after  an  internal  measurement 
and  evaluation  of  the  actual  dynamic  response,  offer  a  more  direct 
method  of  achieving  a  consistent  oontrol  performance.  Among  the 
many  self-adaptive  aircraft  control  systems  whioh  have  been 
proposed  and  investigated  the  Minneapolis-Honeywell  and  General 
Electric  systems  have  probably  undergone  the  furthest  development. 
The  first  of  these  (Ref.  12)  was  first  flown  in  the  X-15  in 
Deoember  1961,  while  a  version  of  the  General  Electrio  system 
(Ref.  6)  is  now  being  installed  in  production  P-111  airoraft. 

Identification 

Inherent  in  the  use  of  self-adaptive  oontrol  systems  is  the 
need  to  identify  the  airoraft  parameters.  An  explicit  identifica¬ 
tion  of  the  aircraft  transfer  funotion  coefficients  involves  & 
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large  amount  of  digital  data  processing  (e.g.  Ref*  22)  on  oontrol 
input  and  output  signals.  The  information  thereby  obtained  is  more 
than  is  required,  unless  the  self-adaptive  system  is  coupled  with 
control  system  parameter  optimization  to  minimize  an  integral 
square  error  criterion  of  performance.  Optimization  in  this  sense, 
simply  to  achieve  acceptable  handling  qualities,  has  not  so  far 
been  found  necessary,  although  in  thrust  or  fuel  management  systems 
such  criteria  are  of  importance  in  minimizing  time  to  climb  to 
altitude,  or  maximizing  range. 

The  success  of  the  Honeywell  and  General  Electric  systems 
undoubtedly  lies  in  the  fact  that  they  oan  be  mechanized  without 
the  need  for  an  airborne  digital  computer.  In  neither  case  is  any 
attempt  made  to  completely  identify  the  aircraft  parameters. 

Instead,  the  assumption  is  made  in  both  systems  that  there  are 

prescribed  regions  in  the  complex  plane  which  define  possible 

locations  for  the  aircraft  zero  and  open  loop  poles.  By  suitable 

choice  of  rate  gyro,  servo/actuator  and  compensator,  the  locus  of 

the  pair  of  complex  actuator  poles  (the  first  to  become  unstable) 

is  forced  to  approach  the  imaginary  axis  in  the  required  direction  / 

for  all  flight  conditions.  The  angle  of  approach  towards  the 

imaginary  axis  is  determined  by  two  different  considerations* 

a.  Honeywell  System.  The  aotuator  poles  approaoh  the  imag¬ 
inary  axis  horizontally,  so  that  the  frequency  of  the  actuator 
mode  (w^)  is  constant  near  the  axis.  A  bandpass  filter  and 
reotif ier  circuit  me  asures  the  amplitude  of  frequency  compon¬ 
ents  at  wQ  and,  in  a  control  loop  containing  a  limiting 
amplifier  and  variable  gain  element,  the  gain  ie  constantly 
adjusted  to  stabilize  a  small  amplitude  limit  cyole.  The 
action  of  the  gain  adjustment  prooess  is  thus  similar  to  that 
in  a  radio  super-regenerative  detector,  and  in  operation  the 
limit  cycle  shows  a  similar  squegging  tendency  (Ref.  1*45)* 

b.  General  Eleotrio  System.  The  aotuator  poles  approaoh  the 
imaginary  axis  with  a  positive  slope  for  increasing  loop  gain, 
so  that  the  loous  passes  through  a  mode  frequency  *ref  with  a 
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damping  factor  of  approximately  0.25*  A  bandpass  filter  cen¬ 
tered  on  is  used  in  conjunction  with  a  frequency  sensor 

to  deteot  on  which  side  of  *re£  'the  actuator  mode  occurs 
whenever  it  is  excited  by  control  Inputs,  and  the  lcop  gain  is 
then  adjusted  to  return  the  mode  to  *Tef 


Scope  of  the  Thesis 

The  self-adaptive  pitch  rate  control  system  which  was  proposed 
and  investigated  in  this  thesis  resembles  the  Honeywell  and  General 
Electric  systems  in  that  complete  parameter  identification  is  not 
made*  only  the  magnitude  of  the  control  surface  effectiveness  is 
determined,  and  this  information  is  used  to  adjust  the  loop  gain  to 
provide  constant  gain  and  phase  margins  for  all  flight  conditions. 

A  distinction  is  that,  whereas  design  of  the  Honeywell  and  General 
Electric  systems  was  carried  out  by  root  locus  methods,  with 
identification  in  the  frequency  domain,  in  this  thesis  the  reverse 
procedure  was  followed!  the  control  loop  design  was  by  frequency 
response  analysis,  and  identification  was  made  in  the  time  domain. 


Figure  1  shows  a  functional  block  diagram  of  the  control 
system,  which  includes  an  airborne  digital  computer  (the  gain 
computer)  and  a  pre-filter  model  reference  to  define  the  desired 
aircraft  response.  The  aircraft  dynamics  were  represented  by  those 
for  the  X-15>  and  the  pre-filter  was  a  first  order  lag  with  a  Time 
constant  of  0.5  sec.  The  way  in  which  the  gain  computer  was  used 
to  obtain  the  required  setting  for  the  variable  gain  element  (Kv) 
is  described  in  the  body  of  the  thesis,  where  results  of  the  system 
simulation  are  presented,  together  with  comments,  conclusions  and 
recommendations. 


Details  of  a  system  simulation  program  for  the  IBM  7094 
computer  are  contained  in  Appendix  A,  and  in  Appendix  B  the  effects 
of  random  wind  gusts  on  the  gain  computer  are  shown.  Approximations 
made  in  deriving  the  aircraft  equations  of  motion  are  discussed  in 
Appendix  C,  while  the  design  cf  a  compensator  for  the  control  loop 
is  considered  in  Appendix  D.  Finally,  difference  equation  approx¬ 
imations  for  the  aircraft  transfer  function  are  given  in  Appendix  E. 
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IX.  Identification  by  Difference  Equations 


The  Aircraft  Transfer  Function 


The  transfer  function  which  relates  aircraft  pitch  rate  (6.) 

B 

to  elevator  deflection  angle  (  S  )  is  derived  in  Appendix  C  as. 


ea(e) 
M  8) 


s  +  l/T 


i  2  2 

8  +2.JW  8  +  wo 


(C.4) 


’a  a 


whor®  the  aircraft  parameter*;  Tft,  and  are  the  short- 
period  time  constant ,  damping  factor  and  natural  frequency,  respec¬ 
tively.  The  parameter  K|  is  the  elevator  effectiveness,  and  from 
the  form  of  equation  (C.4)  it  is  seen  to  be  an  aerodynamic  gain 
fuctor.  In  Appendix  D  it  is  shown  that  a  constant  pitoh  rate 
response  characteristic  may  be  achieved  for  the  control  system  by 
adjusting  only  a  variable  gain  element,  to  compensate  for  changes  in 
the  loop  gain  caused  by  variations  in  Mj .  Thus,  the  self-adaptive 
control  problem  becomes  one  of  identifying  the  elevator  effective¬ 
ness,  and  tracking  changes  in  it  over  a  range  of  values  of  240*1. 


Values  for  the  Airora^t  Parameters 


Figures  2,  3  and  4  illustrate  *’»e  wide  range  of  variation  of 
the  aircraft  parameters  for  thirty-three  Flight  Conditions  (F.C.) 
along  a  typical  flight  path  of  the  X-15»  These  parameters  are 
plotted  against  the  deoibel  value  of  elevator  effectiveness 
relative  to  that  for  F.C.  32,  the  landing  condition.  Specifically, 
for  an  arbitrary  F.C.  n, 

M,  (for  F.C.  n  ) 

-  20  !oc10  (h-  p.c  32) 


0) 


The  usefulness  of  presenting  the  data  in  this  way  is  tha%  as 
•n  approximation,  the  parameter  values  may  be  considered  to  have 
only  one  degree  of  freedom  (Ref.  10 ).  Therefore,  whe»  straight 
lines  are  drawn  through  the  centers  of  the  distributions  ox 
parameter  values  in  Figs.  2,  3  and  4,  the  following  approximate 
relations  are  obtained  if  is  in  dBli 

1/Tft  0.002  exp(0. 152  M*  )  (2) 
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Fig.  4 


Approximation  of  the  Short-Period  *\aturnl  Frequency 
M  •  Function  of  the  Elevator  Effectives .-s* 
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/'+0.50  +  0.03  M,  ?  11  >  M4 

wa  -  V  -0.57  +  0.121  M,  *  ll^Mt^45.5  (3) 

^-40.6  +  Mj  t  M,  >45.5 

Broken  lines  in  Pigs.  2  and  3  are  drawn  for  values  of  l/T&  and 
»w  twice  and  one-half  those  given  by  equation  (2).  Similarly,  in 

Sr 

Fig.  4  "the  broken  lines  are  for  values  of  wa  at  jt20  per  cent  from 
the  values  given  by  equation  (3).  In  the  Figures  it  is  seen  that 
almost  all  of  the  distributions  of  parameter  values  fall  within  the 
pairs  of  broken  lines,  so  that  these  lines  indicate  the  order  of 
approximation  involved  in  estimating  the  other  parameter  values 
from  M|  (in  dBl).  To  test  if  there  was  any  significant  bias  error 
in  equations  (2)  and  (3),  a  calculation  was  made  to  find  the  mean 
errors  involved  in  estimating  1/T  ,  4  w  and  w  from  BL  .  This  was 
done  by  using  the  known  valuer  of  the  parameters  at  each  of  the 
plotted  points.  It  was  found  that  the  mean  error  in  1/T  was  10.9 

8r 

per  cent,  in  ^  w  was  -3.65  per  cent  and  in  w  was  4*62  per  cent. 
Because  of  the  comparatively  much  larger  spread  in  errors  no  attempt 
was  made  to  improve  the  accuracy  of  equations  (2)  and  (3). 

Difference  Equations 

In  Appendix  E  a  method  for  approximating  the  aircraft  pitcrh 
rate  transfer  function  by  a  number  of  finite  difference  equations 
is  presented.  Comparative  tests  with  five  types  of  difference 
equations  showed  (Ref.  Pigs.  E. 3  and  E, 4)  that  1st.  order  central 
differences  gave  root-mean-square  (r.m.s.)  errors  two  orders  of 
magnitude  less  than  the  other  1st.  order  derivative  approximations, 
and  one  order  lees  than  that  of  2nd.  order  backward  differences. 
Therefore,  the  following  oentral  difference  equation  was  chosen  to 
represent  the  aircraft  transfer  function  in  the  gain  computer  (Fig.l), 

T‘*  ■  rrfe 

t-K  T?  M*  1 

*-f~  h\  <B-2> 

ft  J 
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In  equation  (E. 2)  T  is  the  iteration  interval,  and  Y.  ,  is 

s  l+I 

the  calculated  pitch  rate  at  the  instant  when  the  elevator  deflec¬ 
tion  angle  is  Aleo,  the  samples  F>i  and  Yi_1 

represent  the  pitch  rate  and  elevator  deflection  at  the  previous 
iteration,  and  the  one  before  that,  respectively. 


To  find  the  r.m.e.  error  inherent  in  the  use  of  equation  (E. 2) 
the  iterated  values,  Yr  (n  -  1,  2,  3,  . ..),  were  compared  with 
analytically  calculated  values  of  pitch  rate  (^n)»  for  60  elevator 
deflection, 

&e(t)  -  S  [l  +  A3exp(-B3t)  ain(C3t  -  D3)]  (E.7) 

Thus,  equation  (E.7)  represents  the  response  of  a  second  order 
transfer  function  to  a  step  input  of  amplitude  S.  The  difference 
equation  normalized  r.m.B.  error  as  a  function  of  the  number  of 
iterations  was  calculated  as, 

£  (Y  -  P  )2 

n  '  n  n'  _ _  ,  , 

— = — - -  *  100%  (4) 

L.  Pf 


When  the  same  typical  values  of  aircraft  parameters  were  used 
in  both  the  central  difference  equation  (with  T  *  0.01  sec.)  end  in 
the  analytical  solution  for  aircraft  pitch  rate,  the  difference 
equation  r.m.s.  error  was  only  0.01  per  cent  after  200  iterations. 

In  practise,  the  aircraft  pitch  rate  must  be  measured  by  a  rate 
gyro,  which  may  be  expected  to  have  an  error  at  least  this  great. 

It  was  therefore  concluded  that  the  inherent  error  in  the  difference 
equation  (E.2),  when  iterated  at  a  suffiently  fast  rate,  was 
insignificant  compared  with  the  errors  that  would  be  introduced 
by  using  approximations  to  the  aircraft  parameters  in  the  equation. 


Identlf ioation 


Consider  the  self-adaptive  control  system  shown  in  Fig.  1.  At 
first,  when  the  identification  process  is  started,  the  gain  computer 
has  no  knowledge  of  the  values  of  the  aircraft  parameters.  It  does, 
however,  contain  approximate  relations  between  the  values  of  the 
aircraft  parameters,  as  given  by  equations  (2)  and  (3).  Suppose 
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that  there  is  only  one  equation  iterated,  in  the  gain  computer  and, 
since  the  value  of  elevator  effectiveness  is  not  known,  let  be 

*y 

some  arbitrary  but  reasonable  value  used  in  this  y-difference 

equation.  Prom  equations  (2)  and  (3),  approximate  values  for  1/T  , 

ft 

$  w  and  w  may  then  be  computed. 

When  the  subscript  i  is  set  equal  to  zero  in  equation  (E.2)  an 

expression  for  Y,  is  obtained  in  terms  of  Y  ,  E  and  Y  , ,  E  , ,  the 
i  o  0  “1  —l 

initial  conditions,  as  well  as  E^ ,  the  latest  sample  of  elevator 

deflection.  Numerical  values  of  Y  and  Y.  are  available  from  the 

o  i 

rate  gyro  samples  of  aircraft  pitch  rate,  so  that  Yq  ■  Pq  and 


To  ensure  that  the  four  initial  conditions  (denoted  now  by  PQ, 
Eq  and  are  available  in  the  gain  computer  for  the  first 

iteration,  it  is  necessary  to  delay  this  iteration  by  two  sampling 
intervals.  Thereafter,  the  second  and  subsequent  iterations  may  be 
made  as  soon  as  soon  as  each  new  sample  of  elevator  deflection  is 
taken,  without  the  need  for  further  rate  gyro  samples, 

p 

As  each  new  value  Y  is  calculated  the  quantity  (Y  -  P  )  is 

n  '  n  n' 

formed,  to  test  the  accuracy  with  which  the  difference  equation 
response  matches  the  actual  aircraft  response.  Consequently,  after 
m  iterations  the  following  sum-square  error  is  produced, 

S  (m)  -  7  (Y  -  P  )  (5) 

y  '  '  nr/  v/ 

n»l 


By  itself,  S  (m)  is  clearly  not  of  any  use  in  the  identifio- 
ation  process ,  since  it  only  indicates  that  there  was  an  error 
between  the  actual  aircraft  parameters  and  those  used  in  the 
difference  equation;  it  does  not  indicate,  for  example,  whether 
M,  was  too  large  or  too  small.  Furthermore,  even  if  the  correot 

•y 

elevator  effectiveness  had  initially  been  chosen  for  M<  the 

*y 

magnitude  of  S  (a)  would  still  not  be  zero,  because  of  the  very 

if 

rough  approximations  contained  in  equations  (2)  and  (3).  Thus,  at 
least  one  further  difference  equation  (the  x-difference  equation) 
must  also  be  iterated  in  the  gain  computer. 
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Suppose  that,  in  the  x-  and  y-difference  equations,  the  ratio 

Ut  t  M.  is  0.5*  1.0,  and  that  after  every  m  iterations  of  both 

*y 

equations  the  eum-square  errors  S  (m)  and  S  (m)  are  formed.  Suppose 

x  y 

also,  for  example,  that  the  initial  value  arbitrarily  choeen  for 

If.  is  smaller  than  the  true  value  of  aircraft  elevator  effective¬ 
ly 

ness  (Mg). 

If  S  >S  after  the  first  m  iterations,  it  may  be  inferred 
x  y 

that  the  parameters  in  the  y-difference  equation  were  a  better 

match  to  the  true  aircraft  parameters  than  those  in  the  x-differenee 

equation,  and  in  particular  that  <  Mj  <H(.  Therefore,  to  get 

x  y 

closer  to  the  true  aircraft  parameters,  both  lf|  and  are 
increased  by  a  small  fraction  of  their  initial  values.  The  diff¬ 
erence  equations  are  initialized  again,  and  iterated  another  m  times, 
using  further  samples  of  elevator  deflection  as  the  forcing  function. 


The  prooess  is  continually  repeated,  until  it  is  found  that 

S  < S  •  At  this  point  it  is  concluded  that,  in  the  previous  iter- 
*  y 

ation  sequence,  the  parameters  in  the  x-difference  equation  were  a 
better  match  to  the  aircraft  parameters  than  those  in  the  y-diff¬ 
erence  equation,  so  that  M<  <M.<  li.  .  It  is  now  necessary  to 

•x  •  *y 

decrease  the  values  of  M.  and  li,  by  a  small  fraction. 

•x  *y 

Thus,  at  this  stage  the  difference  equations  may  be  "said  to  be 
' looked-on'  to  the  true  elevator  effectiveness.  At  the  end  of  every 

suooeeding  iteration  eequenoe  M.  and  H>  are  incrementally  increa- 

*x  *y 

sed  if  S  >  S  and  decreased  if  S  <  S  to  maintain  the  relationship 
x  y  x  y 

li(  <  H(  «ellj^.  In  this  way  the  value  of  may  be  identified  as 

lying  approximately  mid-way  between  the  values  of  li.  and  M.  •  Also, 

•x  *y 

by  the  use  of  this  type  of  deoision  prooess,  slow  changes  in  the 
magnitude  of  the  aircraft  elevator  effectiveness  will  be  'tracked' 
by  the  difference  equations. 


Error  Correction 

A  disadvantage  in  the  use  of  only  two  difference  equations  for 
identifioation  is  that  the  preoise  value  of  is  not  determined. 
This  problem  may  be  ciroumvented  by  using  an  additional  equation 
(the  z-difference  equation).  In  addition,  the  use  of  three  differ- 
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ence  equations  provides  redundant  information  that  may  be  used  to 

detect  errors  in  the  calculation  of  the  mean-square  errors.  For 

example,  because  of  the  approximations  involved  in  equations  (2) 

and  (3),  it  is  possible  that  sometimes  S  will  be  less  than  S  even 

x  y 

though  M.  may  be  closer  to  M.  than  is  M,  . 

y  *x 

In  the  system  simulation  program  (Ref.  Appendix  A)  normalized 

r.m.s.  errors  were  calculated  for  three  difference  equations  in 

which  If.  <  M,  -c  M.  .  The  following  equation  shows  how  the  r.m.s. 
•x  *y 

error  for  the  y-difference  equation  was  calculated,  and  the  r.m.s. 
errors  for  the  x-  and  z-differenee  equations  were  calculated  in  a 
similar  manner, 

e  m  am 

2 


V"0 


[m  0  /  m 

n«l  /  n-1 


V 


11/2 


(6) 


If  the  r.m.s.  errors  and  the  sum-square  errors  are  each 
arranged  in  ascending  order  after  any  iteration  sequence,  it  may 
he  seen  that  both  sets  will  have  the  subscripts  x,  y  and  z  in  the 
same  relative  positions.  Thus,  either  error  criterion  may  he  used 
to  test  which  difference  equation  response  is  the  best  match  to  the 
actual  aircraft  response.  In  practise,  equation  (5)  is  preferred 
to  equation  (6),  beoause  it  involves  less  computation.  However,  in 
the  system  simulation,  equation  (6)  was  used  because  it  enabled  an 
order-of-magnitude  comparison  to  be  made  between  sets  of  r.m.s. 
errors,  independently  of  the  magnitude  of  the  airoraft  response. 

With  three  difference  equations  iterated  in  the  gain  computer, 
the  following  decision  logic  was  used  to  search  for  and  track 
changes  ir  "ireraft  elevator  effectiveness! 


The  values  of  M 


»x’ 


M,  and  M.  were  all 
’z 


incremented  by  0.6  dB  relative  to  their  previous  values  (a 
factor  of  1,072). 


b. 


Rx<R_2L<R2 


The  values  of  M.  ,  M. 

*x  *' 


and  M.  were  all 
*z 


reduced  by  0.3  dB  relative  to  their  previous  values  (a  factor 
of  1.035). 
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k 


c.  (Ry  <Rx)  AND  (Ry<  R^).  With  this  result,  the  difference 
equations  were  presumed  to  be  lccked-on  to  the  true  elevator 
effectiveness,  and  so  no  change  was  made  in  the  values  of 


I,  ,  m4 

"y 


and  M. 


d.  Other  Combinations.  The  only  other  possible  combinations 
of  the  r.m.s.  errors  are; 

(1)  Ss*iV<Ky>  “*d> 

(2)  Rx<Ss<Hy. 

When  either  of  these  combinations  occurred  the  implication  was 

that  both  Mi  and  Mt  were  closer  to  the  true  Mt  than  was  the 

-  *x  *z  * 

value  of  M.  .  Since  this  is  impossible,  such  results  were 
*y 

ignored,  and  no  changes  were  made  in  the  values  of  U4  ,  U( 

and  M.  .  The  recognition  and  rejection  of  these  combinations 

•a 

of  r.m.s.  errors  constituted  a  partial  check  against  errors 
introduced  in  the  difference  equations  by  the  approximations 
made  in  equations  (2)  and  (3)* 


In  summary,  the  decision  logic  which  was  used  to  change  the 

values  of  M.  ,  M.  and  M,  was  designed  with  the  object  of  making 
•x  *y  *2, 

the  r.m.s.  error  associated  with  the  y-difference  equation  the 
smallest  of  the  three.  Consequently,  the  identified  value  of 

elevator  effectiveness  was  taken  to  be  equal  to  the  value  of  M.  . 

•y 

The  variable  gain  element  (Ky)  in  the  control  loop  was  then  changed 

inversely  to  changes  in  ,  as  given  by  equation  (D. 6).  Thus, 

when  M,  was  increased  by  0.6  dB  the  value  of  K_  was  decreased  by 
•y  ▼' 

0.6  dB,  and  when  Mt  was  decreased  by  0,3  dB  the  value  of  K  was 

y 

increased  by  0.3  dB.  The  step  decreases  in  Ky  were  made  larger 
than  the  increases  to  ensure  stability  of  the  control  loop  in 
flight  conditions  where  the  aircraft  elevator  effectiveness  was 
increasing. 


Length  of  Iteration  Sequence 

A  further  limitation  introduced  by  the  use  of  equations  (2) 
and  (3)  is  that  if  the  iteration  sequenoes  are  too  long  the  r.m.s. 
errors  beoome  so  large  that  identification  of  Mj  is  Impossible. 
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Before  a  simulation  program  for  the  complete  closed  loop  control 
system  was  written,  an  investigation  was  therefore  made  to  find  a 
suitable  iteration  length.  In  a  similar  manner  to  that  previously 
described  (Ref.  page  11)  for  finding  the  error  inherent  in  the  use 
of  difference  equations,  r.m.s.  errors  were  computed  for  the  x-,  y- 
and  7, -difference  equations  as  a  function  of  the  number  of  iterations. 

A  value  of  0.01  sec.  was  used  for  x  ,  and  the  ratios  of  If.  i  If.  t  If. 

s’  *x  *y 

were  0.5  *  1.0  t  2.0. 

Analytical  values  for  the  aircraft  response  to  the  elevator 
deflection  function  given  by  equation  (E. 7)  were  calculated  for 
ten  flight  conditions.  For  the  same  elevator  deflection  function, 
r.m.s.  errors  for  the  x-,  y-  and  z-differenee  equations  were 
calculated  for  values  of  If.  removed  from  the  true  value  of  elevator 

*y 

effectiveness  by  up  to  +30  dB  and  -36  dB.  To  provide  a  represen¬ 
tative  sample  of  results,  twenty-one  computer  runs  were  made  in 

which  the  combinations  of  flight  condition  and  value  of  Mg  relative 

y 

M(  were  chosen  at  random. 

For  values  of  m  equal  to  10,  20,  40,  60,  80  and  100  the 

decision  rules  for  incrementing,  reducing  or  leaving  unchanged  the 

values  of  M.  ,  Mg  and  Mg  were  applied  to  the  r.m.s.  errors  R  (o), 

Ry(m)  and  Rg(m).  A  correct  decision  was  assumed  to  have  been  made 

if  the  r.m.s.  errors  indicated  a  need  for  a  change  in  Mi  which 

,  y 

would  have  made  it  more  nearly  equal  to  Mg  (or  left  Mg  unchanged 
if  it  waa  already  equal  to  Mg). 

The  reeulte  of  this  investigation  are  presented  in  Fig.  5» 
where  it  ie  shown  that  the  r.m.s.  error  oriterion  gave  21  oorreot 
deoisions  out  of  21  up  to  an  iteration  length  of  0.4  seo.  (i.e.  up 
to  40  iterations).  With  longer  iteration  times  the  number  of 
oorreot  decisions  was  reduoed,  until  after  one  stcond  only  14 
oorreot  deoisions  were  made. 

In  addition  to  r.m.s.  errors,  the  correlation  ooeffioients 
were  oaloulated  between  the  difference  equation  responses  and  the 
airoraft  response.  The  following  equation  defines  the  correlation 
ooeffioient  between  the  y-difference  equation  response  and  the 
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^correct 
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(tut  of  ll) 
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Iteration  length 


Fig.  5 


Performance  of  the  Decision  Logiv-  Whioh  wag  Used  to  bewrch 
for  and  Traok  Changes  in  Elevator  Effect ivenees 
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aircraft  response,  and  the  correlation  coefficients  for  the  x -  and 
«-difference  equations  were  calculated  in  a  similar  manner, 

|l/2 


C  (m)  » 

yv 


m  >  r  m  „  m 

r  tr  p ) / 1  E Y  If 

n«l  n  n  /  ln-1  n  n-1  nJ 


(7) 


If,  in  equation  (7)  all  the  difference  equation  samples,  Yn, 
are  multiplied  by  a  constant  factor,  the  value  of  C  (m)  is  unohanged. 
The  values  of  C^.,  Cy  and  are  not,  therefore,  directly  affected  by 
the  different  values  of  elevator  effectiveness  used  in  the  x-,  y- 
and  z-difference  equations.  Consequently,  the  correlation  ooeffio- 
ients  depend  only  on  the  fidelity  with  which  the  wave-shapes  of  the 
difference  equation  responses  match  the  aircraft  response.  The 
magnitudes  of  the  correlation  coefficients  are  particularly  sensitive 
to  mismatches  in  the  times  at  which  axis  crossings  occur.  It  follows 
from  these  considerations  that  the  degree  of  correlation  depends 
only  on  the  accuracy  of  the  approximations  to  the  aircraft  param¬ 
eters  1/T  ,  jr w  and  w  used  in  the  difference  equations  (for  a 
given  elevator  deflection  function). 

In  Pig.  5»  a  correot  decision  is  plotted  for  the  oaseB  where 
the  difference  equation  which  best  correlated  with  the  aircraft 
response  also  used  an  elevator  effeotivenese  nearest  in  value  to  • 
Beoause  there  were  three  difference  equations,  a  correot  decision 
had  a  probability  of  l/3rd.  of  ooourring  by  ohanoe.  Up  to  the  first 
0.1  sec.  iteration  of  the  difference  equations,  their  correlations 
with  the  airoraft  response  were  good  enough  to  give  20  oorreot 
decisions  out  of  21.  Thereafter,  when  the  difference  equation 
responses  started  to  diverge  from  the  airoraft  response,  the  number 
of  oorreot  deoisions  was  worse  than  chance  expectation.  Towards  the 
end  of  the  run,  however,  the  decision  by  correlation  coefficients 
began  to  improve  slightly,  as  the  values  of  the  ooeffiolents 
stabilized  with  the  increased  sample  else. 

Based  on  the  results  of  this  investigation,  an  iteration 
length  of  0.1  sec.  was  seleoted  for  the  identification  interval. 

It  was  also  found  that  when  ■  M(  the  value  of  Be  was  consistently 
greater  than  so  that  the  ratios  M|x  t  M4y  i  were  ohanged  from 
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0.  5  t  l.Ot  2.0  to  the  ratios  0.5*  1.0*  1.5  for  the  closed  loop 
system  simulation  program. 

Closed  Loop  Test  of  ths  Identification  Process 

A  simulation  program  (Ref.  Appendix  A)  was  written  for  the 
IBM  7094  computer  to  investigate  the  usefulness  and  accuracy  of  the 
identification  process.  Thij  program  simulated  the  complete  self- 
adaptive  pitch  rate  control  system  shown  in  Pig.  1.  It  included 
the  X-15  dynamics  for  four  widely  different  flight  conditions 
(Ref.  Appendix  C,  Table  II),  all  the  control  system  components, 
and  the  airborne  digital  computer. 

Results  of  +he  simulation  for  pilot  command  inputs  are  given 
in  Chapter  IV;  results  for  random  wind  gusts  are  given  in  Appendix  B, 
together  with  details  of  the  technique  used  to  generate  band-limited 
noise  samples  to  represent  wind  gusts.  The  objectives  of  the 
investigation  were: 

a.  To  demonstrate  that  the  high-gain  control  loop  was  capable 
of  forcing  the  aircraft  pitch  rate  to  accurately  follow  the 
output  of  the  pre-filter  model  reference. 

b.  To  show  that  by  adjusting  a  variable  gain  element  in  the 
loop  to  compensate  for  changes  in  elevator  effectiveness  (Mj ) 
a  consistently  accurate  control  performance  could  be  achieved 
for  all  flight  conditions. 

c.  To  test  the  speed  and  accuracy  with  which  the  airborne 
computer  (the  gain  computer)  was  able  to  search  for  and 
identify  the  airoraft  elevator  effectiveness,  and  set  the 
control  loop  gain  accordingly  (Ref.  equation  (D.6)  ). 

d.  To  confirm  that  the  identification  and  self-adaptive 
process  was  feasible,  without  any  need  for  external  air  data 
information. 

e.  To  investigate  the  effect  of  severe  wind  gusts  on  the 
performance  of  the  eelf-adaptive  system. 
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III.  Simplified  Analysis  of  the  Identif ioation  Process 
The  Simplified  Control  Loop 

Pilot  pitch  rate  commands  are  passed  through  a  pre-filter, 
which  is  a  first-ordsr  lag  with  a  break  frequency  of  2  rad. /sec. 

Por  this  analysis  it  is  therefore  assumed  that,  for  pilot  commands, 
the  signal  components  in  the  control  loop  are  mainly  composed  of 
low  frequency  terms.  Thus,  as  a  first  approximation,  all  transfer 
functions  are  replaced  by  their  low  frequency  gain  factors,  as 
shown  in  Fig.  6. 


It  is  assumed  that  the  ratios  of  the  difference  equation  gain 

factors,  K  «  K  i  K  are  given  by  0.5«  1.0 »  1.5*  Since  a  positive 
X  jr  * 

pilot  command  Input  oauses  a  negative  elevator  deflection,  which 
in  turn  results  in  a  positive  aircraft  pitoh  rate,  all  the  gain 
faotora,  K  ,  K  ,  K  ,  K  and  K  ,  are  positive.  Term#  with  a  auV 

O  o  X  j  ■ 

soript  n  are  sampled  values  of  signals  In  the  oontrol  system  at  the 

difference  equation  iteration  times  (every  seo.). 

0 
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Prom  equation  (5), 


s 

z 

m 

e<vv2 

n 

(8a) 

sy 

* 

E<W2 

n 

(8b) 

Sx 

* 

e<vp„)2 

n 

(8c) 

But,  Zw  -  -K  E  , 
n  an 

y 

n 

■  -K  E  and  X  ■  -  K  E  , 
y  n  n  x  n 

so  that. 

S 

z 

- 

r  (k  e  +  p  'i2 

'  z  n  n 

n 

(9a) 

S 

y 

m 

pyW2 

(9b) 

s 

X 

- 

n 

(9c) 

Case  1  -  Command  Inputs  Only 


In  this  case,  «  0 

for  all  n,  and 

Pn 

•  -K  E 
a  n 

(10) 

Substituting  in  equations  (9)  for  Pn, 

S  « 

z 

<VV2E** 

n 

(Ha) 

sy  ■ 

(Ky“Ke>2  EEn 
n 

(Hb) 

« 

(K  -K  )2  Fe2 

x  a  “  n 
n 

(He) 

When  the  difference  equations  are  locked-on  to  the  true 

value 

elevator  effectiveness, 

K  ■  K  ,  and  also, 
y  a  ’ 

■  0-5Kt 

end 

'  1-5K. 

Substituting  in  equations  (11)  for  K^, 

Ky 

and  R  , 
z 

Sz  * 

«•»<  E< 

n 

(12a) 

sy  " 

0 

(12b) 

Sx  ■ 

°-»*2E*2 

(12c) 

Thus,  when  the  difference  equations  are  locked-on  and  the 
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gain  factors  of  these  equations  are  in  the  ratios  0.5 *  1.0*  1.5* 

the  sum-square  errors  (and  the  r.m.s.  errors)  are  symmetrical  with 

respect  to  the  y-difference  equation  error,  which  ideally  is  zero. 

However,  the  behavior  of  the  errors  is  different  when  K  >K  than 

y  a 

when  K  <K  .  For  example,  if  K  is  too  large,  say  K  -  10K  ,  then 
ya  y  y  a 

substitution  of  K  and  K z  into  equations  (11)  gives  ratios  for 
the  r.m.s.  errors, 


R  *  R  *  R 

x  y  z 


equal  to  4*9*14 


But,  if  K  is  too  small,  say  K  *=  0.1  K  ,  then  the  ratios  of 

y  y  a 


the  r.m.s.  errors  are, 


R  *  R  t  R 
x  y  z 


equal  to  19  *  1$  *  17 


Hence,  when  the  y-difference  equation  gain  is  greater  than  the 

aircraft  gain  the  r.m.s.  errors  are  more  widely  separated  than  when 

the  y-difference  equation  gain  is  less  than  that  of  the  aircraft. 

The  decision  logic  (Ref.  page  14)  is  consequently  more  likely  to 

detect  correctly  a  requirement  to  reduce  Mj  (and  increase  Ky) 

than  to  increase  M,  (and  reduoe  K  ). 

<y  v 

After  this  behavior  of  the  r.m.s.  errors  was  substantiated  in 
the  closed  loop  simulation  of  the  system,  advantage  was  taken  of  it 
to  ensure  that  Kv  would  not  unnecessarily  be  increased,  thereby 
reducing  the  control  system  stability.  To  do  this  the  decision 

logic  to  reduoe  M.  was  changed  from  the  requirement  that 

•y 

R  c  R  «*R  to  the  stronger  requirement  that  (R  -<R  )AND(3R  <R  ). 
x  y  z  y  z  x  z 

The  results  of  this  change  are  shown  in  the  next  chapter. 


Case  2  -  Disturbance  Inputs  Only 

Although  the  primary  effect  of  vertical  wind  gusts  is  to  cause 
a  ohange  in  the  aircraft  angle  of  attaok  (Ref.  Appendix  B),  the 
only  interaction  with  the  control  system  is  through  the  resulting 
aircraft  pitoh  rate.  Therefore,  in  this  analysis  the  wind  gusts 
ars  considered  in  terms  of  a  disturbance  pitoh  rate,  samples  of 
whioh  are  represented  in  Fig.  6  by  Gr. 

For  no  pilot  oommands,  ■  0  for  all  n. 

Henoe,  ■  +K  P_  (13) 

non 
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As  in  the  previous  case,  when  the  difference  equations  are 
locked-on,  K  -  0.5K  .  K  ■  K  and  K  -  1.5K  . 

JL  B  j  91  Z  a 

Substituting  in  equations  (?)  for  E  ,  K  ,  K  and  K  ,  the  sum- 

ft  x  y  z 

square  errors  now  become, 

S.  *  <1-5fa,to-tl)2Ern 

n 

sy  •  <KaKo  +  1>2IPn  (^b) 

ft 

s*  ■  <°-5K.!C0  +  l)2£pp  U4C) 

n 

Sinoe  K  and  K  are  positive,  S  <  S  <:  S  .  With  this  result. 

«  o  x  y  z 

the  decision  logic  provides  for  the  gain  faotor  of  the  y-differenoe 

equation  to  be  reduced.  This  analysis  therefore  shows  that  low 

frequency  pitch  rate  disturbances  tend  to  increase  the  value  of  K  . 
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IV.  Performance  of  the  Gain  Computer  in  the  System  Simulation 

Figures  7  through  20  illustrate  the  performance  of  the  self- 
adaptive  control  system  for  pilot  command  inputs.  These  results 
were  obtained  by  means  of  a  simulation  program  which  was  run  on 
an  IBM  7094  computer.  Because  of  the  large  number  of  integrations 
required  to  simulate  the  system  transfer  functions,  the  program 
speed  was  approximately  l/30th.  of  real  time.  The  results  of  the 
computations  were  stored  on  magnetic  tape,  and  then  subsequently 
plotted  on  a  Benson-Lehner  digital  plotter.  Details  of  the  program 
are  contained  in  Appendix  A. 

The  quantities  plotted  on  the  graphs  are: 

a.  Pilot  command  input,  ©  (deg. /sec.). 

c 

b.  Output  of  the  pre-filter  model  reference,  0m  (deg. /sec.). 

c.  Aircraft  pitch  rate,  ©  (deg. /sec.). 

ft 

d.  Elevator  deflection  angle  $e  (deg.). 

e.  Gain  error  in  the  control  loop  relative  to  the  design 

value  (Ref.  equations  (D.6)  and  (D.7)  )  for  the  particular 

Flight  Condition  (F.C.),  dB. 

For  each  F.C.  the  loop  gain  was  initially  offset  from  its 
design  value  in  order  to  test  the  speed  and  accuracy  with  which  the 
gain  computer  returned  it  to  the  design  value.  The  amount  of 
initial  offset  was  in  some  cases  restricted.  For  example,  in 
F.C.  28  the  design  value  of  loop  gain  was  achieved  when  the 
variable  gain  element  was  at  its  minimum  setting,  so  that  the  loop 
gain  could  never  he  too  low  for  this  F.C.  Since  the  oontrol  loop 
had  a  design  value  of  gain  margin  of  approximately  12  dB  for  all 
F.C.,  it  was  not  practicable  to  set  in  a  positive  gain  error  higher 
than  this.  Thus,  in  Fig.  8  the  control  loop  is  seen  to  be  just  on 
the  verge  of  instability  until  the  gain  error  is  reduced  below  12  dB. 

In  all  cases  when  the  loop  gair  setting  was  not  gro«»?ly  in 
error  the  aircraft  response  followed  the  model  response  with 
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negligible  error.  By  comparing  the-  elevator  motion  for  F.C.  28 
in  Fig.  7  vith  that  for  F.C.  32  in  Fig.  12,  the  wide  variation  in 
elevator  effectiveness  between  the  F.C.  is  apparent. 

The  decision  logic  which  wee  used  in  the  gain  computer  to 
obtain  Figures  7  through  12  was  that  which  is  described  on  page  14. 
With  the  exception  of  F.C.  7  (Fig.  10),  the  gain  changing  was 
erratic,  although  the  final  gain  errors  were  small?  4  dB  for 
F.C.  28,  and  less  than  2  dB  for  the  other  flight  conditions. 

When  a  small  amplitude  dither  signal  wae  injected  into  the 
control  loop  to  continuously  excite  the  system  response,  the 
performance  of  the  gain  computer  was  considerably  improved.  A 
comparison  of  Fig.  13  with  Fig.  7  shows  that  with  the  dither 
signal  there  was  less  tendency  for  the  gain  to  stay  fixed  at  too 
high  a  value.  This  is  particularly  noticable  in  the  first  two 
seconds  of  each  run,  when  there  was  no  pilot  input.  In  Fig.  14, 
the  effect  of  modifying  the  decision  logic  (Ref.  page  22)  to 
prevent  unnecessary  increases  in  gain  is  shown.  Thus,  both  of 
these  changes  improved  the  performance  of  the  gain  computer  for 
F.C.  28. 

With  theee  two  changes,  repeats  of  the  computer  runs  for  the 
other  three  F.C.  were  then  made.  A  comparison  of  Figs.  8  to  12 
with  Figs.  15  to  19,  respectively,  shows  that  in  each  case  the 
gain  errors  were  reduced,  and  gain  adaption  took  place  more  rapidly. 

The  dither  signal  whioh  mas  used  in  these  tests  had  a  frequency 
of  30  rad. /sec.  and  was  inserted  into  the  control  loop  at  the  input 
summing  junotion.  To  prevent  the  olevator  motion  from  becoming 
excessive  in  conditions  of  low  elevator  effectiveness,  the  amplitude 
of  the  signal  was  progressively  reduoed  as  the  identified  value  of 
elevator  effectiveness  became  smaller.  For  F.C.  28,  the  condition 
of  greatest  elevator  effectiveness,  the  amplitude  was  0.01  deg./seo. , 
whereas  for  F.C.  32,  the  oondition  of  least  elevator  effectiveness, 
the  amplitude  was  reduoed  to  l/5th,  of  this  value.  Fig.  20  shows 
the  exosssive  elevator  motion  whioh  ooourred  when  a  variable  amplit¬ 
ude  dither  signal  was  not  used. 
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Fig.  7 

Flight  Condition  28  -  Loop  Gain  Initially  12  dB  Too  High 
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Pig.  10 

Plight  Condition  7  -  Loop  Gain  Initially  16  dB  Too  Low 
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Fig.  20 

Flight  Condition  32  -  Showing  Excessive  Elevator  Motion 
a  Constant  Amplitude  Dither  Signal  was  Used 
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V.  Cone  1  as  ions  and  Recommendations 


Conclusions 

•1.  Identification  Accuracy.  Difference  equations,  iterated  in 
an  airborne  digital  computer,  may  be  used  to  identify  the  value  of 
elevator  effectiveness  within  approximately  ±2  dB,  provided  that 
wind  gust  disturbances  are  not  present  and  that  accurate  measure¬ 
ments  of  elevator  deflection  and  aircraft  pitch  rate  are  available. 

2.  Loop  Gain  Adjustment.  By  adjusting  a  variable  gain  element 
in  the  forward  loop  inversely  to  changes  in  the  identified  value 
of  elevator  effectiveness  the  gain  crossover  frequency  may  be  held 
within  the  range  of  27  to  37  rad. /see.  This  frequency  range  is 
only  0.14  of  a  decade.  In  comparison,  the  gain  crossover  frequency 
shifts  over  approximately  1.7  decades  for  a  fixed  gain  system, 
which  greatly  increases  the  problem  of  stabilizing  the  control 
loop.  The  figure  of  1.7  decades  is  calculated  by  assuming  that, 
for  stability,  an  average  gain  slope  of  28  dB/decade  is  required  in 
the  region  of  the  gain  crossover  frequency,  which  shifts  as  the 
elevator  effectiveness  varies  over  47*6  dB  (240il). 

3.  Bandwidth.  The  design  open  loop  bandwidth  (defined  by  the 
gain  crossover  frequency)  for  the  self-adaptive  system  was  set  at 
31.3  rad. /sec.  At  this  value,  a  fast  oontrol  response  was  aohieved, 
and  yet  the  bandwidth  was  not  so  high  that  body  bending  modes 
would  be  significant.  If  a  fixed  gain  system  had  been  used,  with 
this  bandwidth  for  Plight  Condition  32  (llj  ■  -0.219)»  then  the 
bandwidth  for  Flight  Condition  28  (Mt  «  -52.9)  at  a  frequency  of 
1.7  decades  higher  would  have  beon  470  rad./seo.  Although  this 

is  an  approximate  theoretical  figure,  it  is  dose  to  the  value  of 
520  rad. /sec.  which  was  found  to  be  necessary  by  Beals  and  Sellings 
(Ref.  1«94).  In  practise,  of  oourse,  such  large  bandwidths  art 
quite  impossible  to  realize  beoause  of  saturation  effects  in 
physioal  components  in  the  oontrol  loop,  noise  piok-up  and  the 
deetabilizing  action  of  bigh  frequency  body  bending  modes.  For  the 
open  loop  bandwidth  of  31.5  rad./seo.  the  closed  loop  bandwidth* 
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measured  from  the  output  of  the  pre-filter  model  reference  to  the 
aircraft  pitch  rate,  was  6  dB  down  at  a  frequency  of  22  rad* /sec. 

4«  Response  Speed.  The  control  system  response  speed  and 
consistency  between  the  extreme  ranges  of  elevator  effectiveness 
is  demonstrated  by  the  following  example.  For  a  step  pilot 
command  input  the  output  of  the  pre— filter  model  reference  reached 
90  per  cent  of  its  steady-state  value  after  1.15  sec.  With  this 
i  command,  the  aircraft  pitch  rate  reached  90  per  cent  of  the  steady- 

state  value  after  1.26  sec.  for  Flight  Condition  28,  and  after 
1.14  sec.  (implying  a  slight  overshoot)  for  Flight  Condition  32. 

5»  Dither  Signal.  The  injeotion  of  a  small  amplitude  dither 
signal  into  the  control  loop  to  continuously  excite  an  aircraft 
response  greatly  improved  the  speed  and  accuracy  of  the  identifica¬ 
tion  process.  The  frequency  of  the  dither  signal  was  30  rad. /sec. , 
and  its  maximum  amplitude  was  0.01  deg./sec.  To  prevent  excessive 
elevator  motion  under  conditions  of  low  elevator  effectiveness,  the 
amplitude  was  reduoed  as  a  function  of  the  identified  elevator 
effectiveness,  so  that  for  Flight  Condition  32  the  amplitude  was 
only  1/5  of  the  maximum  value.  In  a  practical  implimentation  of 
the  identification  process  cirouit  noise  and  small  disturbances 
to  the  aircraft  flight  path  would  probably  be  sufficient  to  make 
identification  possible  without  a  dither  signal.  This  was  found 
to  be  true  for  the  General  Eleotrio  system  (Ref.  6*89). 

6*  liSi  Gu8t8«  presence  of  wind  gusts  degraded  the  accuracy 

with  which  the  elevator  effectiveness  was  identified.  For  severe 
random  wind  gusts  with  an  r.m.s.  magnitude  of  20  ft. /sec.  their 
effeot  on  the  oontrol  system  was  to  oause  the  loop  gain  to  decrease 
by  4  to  6  dB  below  its  design  value.  Because  of  the  tolerance  of 

* 

the  high  gain  feedback  loop  to  gain  variations  of  this  size,  there 
was  no  adverse  effect  on  the  control  response  characteristics. 

7.  Differenoe  Equations.  To  provide  the  greatest  acouraoy  when 
transfer  functions  are  represented  by  finite  difference  equations, 
derivative  approximations  should  be  made  by  means  of  central 
difference  formulas.  For  a  given  test  funotion,  central  differences 
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were  found  to  give  r.m.s.  errors  of  only  0.01  per  cent,  compared 
with  errors  of  approximately  4  Per  cent  for  both  backward  and 
forward  differences. 

Recommendations. 

1.  A  low  frequency  analysis  of  the  behavior  of  the  decision  logic 
indicated  that  wind  gust  disturbances  acting  on  the  aircraft  should 
cause  the  loop  gain  to  be  increased  beyond  its  design  value. 
However,  in  the  system  simulation  the  gain  was  found  to  decrease. 
Because  of  lack  of  time,  this  phenomenon  was  not  extensively 
investigated,  although  it  appears  likely  to  be  due  to  failure  of 
the  difference  equations  to  respond  accurately  to  high  frequency 
components  in  the  elevator  deflection  motion.  It  is  therefore 
recommended  that  the  effect  of  wind  gusts  on  the  performance  of 
the  gain  computer  should  be  the  subject  of  further  investigation. 

2.  The  central  difference  equation  approximation  to  the  aircraft 
transfer  function  was  found  to  be  very  accurate  (0.01  per  cent 
r.m.s.  error)  when  the  exact  aircraft  parameters  were  used  in  it 
for  test  purposes.  It  is  therefore  recommended  that  a  study  be 
made  to  investigate  the  feasibility  of  writing  a  genexal  purpose 
computer  program  for  the  simulation  of  closed  loop  control  systems 
on  the  AFIT  IBM  1620  machine.  Such  a  program  would  use  a  oentral 
difference  equation  to  represent  each  transfer  function  in  the 
control  loop. 
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Appendix  A 

The  System  Simulation  Program 


■Input  Data 

Input  data  for  the  program  is  contained  on  three  cards.  The 
parameters  which  ar®  required  to  he  specified  on  these  cards  are 
listed  below. 


Card  No.  1 

WVEL  A  velocity  (ff„  )  as  a  starting  value  for  the 

wind  gust  samples  defined  in  equation  (B.8). 

UAC  Aircraft  velocity  (U)j  specified  in  Appendix  G, 

Table  II. 

TS  Integration  step  interval  (t  ),  and  also  wind 

8 

gust  sample  interval. 


Card  No.  2 

GP 

GV 

STOP 

K 


TEMPG 
Card  No.  3 
MQ 

MALF 

MDEL 

LALF 

LDEL 


Control  loop  fixed  gain  factor  (Kf). 

Initial  value  for  the  variable  gain  factor  (K  ). 

V 

Length  of  simulation  run  (sec.)  required. 

A  switch.  If  K  ■  1.0  the  program  will  run 
normally,  but  if  K  »  0.0  then  Kv  will  remain 
fixed  at  its  initial  value. 

A  dummy  variable,  which  may  be  ignored. 


The  aircraft  stability  derivatives  Mq, 

Mj ,  Lft,  and  Lj ,  respectively.  Numerical 
values  are  contained  in  Appendix  C,  Table  I. 


Program  Organization 

A  list  of  the  program  statements  is  given  at  the  end  of  thic 
Appendix.  It  will  be  noticed  that  the  complete  pregram  is  divided 
into  two  sections j  the  first  is  written  in  the  illNIC  rrogrammlng 
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language,  and  the  second  in  FORTRAN.  In  general  terras,  the  two 
sections  of  the  program  perform  the  following  operations t 

MIMIC 

a*  Simulation  of  the  aircraft  longitudinal  equations  of 
motion;  equations  (C.l)  and  (C.2). 

b.  Simulation  of  the  components  comprising  the  control 


system; 

(1) 

Pre-filter  model  reference. 

(2) 

Fixed  and  variable  gain  elements 

(3) 

Compensator. 

(4) 

Servo  and  aetuator. 

(5) 

Rate  gyro. 

o.  Generation  of  inputs  to  the  control  system; 

(1)  Pilot  pitch  rate  oommands. 

(2)  Dither  signal  to  stimul&Le  the  adaptive  control 
mechanism  in  the  absence  of  pilot  oommands. 

(3)  Random  wind  gusts.  (Ref.  Appendix  B) 

d.  Transmit xot,  of  data  to  the  airborne  digital  computer 
(the  gain  computer,  simulated  in  the  FORTRAN  section); 

(1)  Samples  of  elevator  deflection  angle  and  aircraft 
pitch  rate  (as  measured  by  the  rate  gyro)  every  0.01  sec. 

(2)  The  present  value  of  the  variable  gain  setting 
every  0.1  sec. 

e.  Output  of  pilot  oomraande  and  aircraft  responses,  etc., 
for  storags  on  magnetio  tape.  The  tape  could  then  be  later 
transcribed  on  the  Beneon-Lehner  digital  plotter. 

FORTRAN 

s.  Simulation  of  the  airborne  digital  computer; 

(1)  Iteration  of  the  three  difference  equations,  and 
calculation  of  their  roet-mean-equare  errors., 

(2)  Decision  whether  or  not  to  change  tho  loop  gain. 

b.  Transmission  of  output  data  to  the  line  printer; 

(l)  Errors  in  the  difference  equation  solutions. 
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(2)  Control  loop  gain  error,  calculated  from  equation 
(D.7).  It  should  he  noted  that  this  calculation  of  gain 
error  (£Kv)  is  readily  made  during  simulation,  where  the 
true  value  of  elevator  effectiveness  for  the  Flight 
Condition  which  is  being  aimulated  is  known  precisely. 

Program  Execution  Sequence 

Through  the  use  of  the  MIMIC  processor  (Ref.  17),  which  is 
available  on  the  IBM  7094  computer,  differential  equations  may  be 
set  up  and  solved  in  a  way  which  closely  resembles  that  used  in 
analog  computation.  Consider,  for  example,  equation  (C.l)t 

0a  “  uq  e>a  *  Ma  a  +  M,  (c.l) 

In  the  MIMIC  program  language  this  equation  is  written, 

DQAC  •  MQ*QAC+MALP*ALF+MDEL*£L 

where,  DQAC  is  the  program  symbol  for  ©  ,  QAC  is  the 

a 

symbol  for  ft  ,  and  the  other  terms  are  self-descriptive.  Then, 

C* 

to  get  &  from  $  ,  the  variable  DQAC  is  simply  used  as  the  input 
to  a  MIMIC  integrator,  with  zero  initial  conditions i 

QAC  »  INT ( DQAC >0* ) 

Unlike  an  analog  computer,  however,  the  digital  computer  can 
only  perform  one  integration  at  a  time,  so  that  the  integrator 
outputs  must  be  updated  in  sequence ?  this  is  organized  automatic¬ 
ally  by  the  MIMIC  processor.  Normally,  the  processor  also  adjusts 
the  integration  step  interval  so  that  the  error  at  the  output  of 
each  integrator  is  held  to  within  5  parts  in  10^.  When  the  random 
wind  gust  inputs  were  used  this  latter  facility  was  not  available, 
and  instead,  a  fixed  integration  step  of  0.0005  sec.  nas  used  as  a 
compromise  between  accuracy  and  speed  of  computation.  This  step 
size  dictated  the  sampling  period  (t  )  for  the  wind  gusts, 

3 

To  provide  computation  time  for  simulation  of  the  airborne 
computer,  the  following  time-sharing  routine  was  established  in  the 
IBM  7094  computer! 
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a.  The  initial  value  for  the  setting  of  the  variable  gain 
element  (Kv),  which  was  input  to  the  program  on  a  data  card, 
wan  used  in  the  MIMIC  simulation  of  the  control  loop  during 
the  first  0.1  sec. 


b.  Every  0.01  sec.  the  computed  values  of  elevator  deflection 
and  aircraft  pitch  rate  were  stored  and  transfered  to  the 
FORTRAN  section  of  the  program. 


c.  After  the  end  of  the  first  0.1  sec.  of  MIMIC  computation 
(and  at  subsequent  0.1  sec.  intervals)  all  the  outputs  of  the 
MIMIC  integrators  were  'frozen1,  and  computation  was  transfersd 
to  the  FORTRAN  section. 


d. 


Using  the  known  value  of  Kv,  the  best  estimate  of  computed 


elevator  effectiveness  (Mr  )  was  established  from  the  relation- 

•y 

ship  between  these  two  variables  defined  by  equation  (D.6). 


e.  Coefficients  for  the  three  difference  equations  were 

calculated  from  M.  ,  as  described  in  Appendix  E. 

*y 

f.  The  difference  equations  were  initialized  and  then  iterated, 
using  the  stored  values  of  elevator  deflection  as  inputs. 

g.  A  calculation  of  the  root-mean-equare  errors  of  the  diff¬ 
erence  equations  was  made  by  comparing  the  iterated  solutions 
with  the  stored  values  of  aircraft  pitch  rate. 


h.  Based  on  the  relative  magnitudes  of  the  p.m.s.  errors 
(R^,  Ry,  and  R^),  *  decision  was  made  whether  or  not  to  change 
the  loop  gain i 

(1)  R<R<R  |  K  increased  by  0.3  dB. 

n/  x  y  z  v 

(2)  R  >  R  >  R  (  K  decreased  by  0.6  dB. 

'  X  y  z  v 

(3)  Rx>  Ry  and  R^  >  Ry  |  Kv  unchanged. 

i .  MIMIC  simulation  for  a  further  0.1  seo.  was  continued, 

with  the  new  value  of  Kv,  and  the  routine  was  repeated  until 
the  required  program  running  time  had  been  completed.  It  may 
be  noted  that  equation  (D.6)  was  always  used  to  relate  Kyr  and 

M.  when  the  commutation  was  transfared  from  the  MIMIC  seotion 

•y 

to  the  FORTRAN  section,  or  vice-versa. 
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Output  Data 

A  short  sample  of  printed  output  data  from  the  IBM  7094  com¬ 
puter  is  given  in  Fig.  A.l.  These  results  were  obtained  by  the 
application  of  a  0.5  deg. /sec.  step  command  input  to  the  oontrol 
system,  with  the  aircraft  in  Flight  Condition  24,  and  with  an 
initial  gain  error  (AKv)  of  6  dB. 

Printed  aoross  the  page  are  values  of  the  aircraft  pitoh  rate 
(from  the  rate  gyro)  at  0.01  sec.  intervals,  and  underneath  them 
are  the  corresponding  errors  in  the  difference  equation  solutions. 

It  can  be  seen  that  the  difference  equations  were  restarted  every 
0.1  sec.,  using  two  samples  of  the  aircraft  pitoh  rate  to  set 
in  the  initial  conditions.  The  extreme  right  hand  column  contains 
the  final  values  of  the  difference  equation  errors  at  the  end  of 
each  0.1  second  computation  cycle  (i.e.  after  eight  iterations). 

Also  in  this  column  is  the  gain  error,  which  decreased  from  6  dB  to 
3  dB  over  the  period  of  0.7  sec.  To  the  left  of  this  column  are 
three  offset  columns  containing  the  root-mean-square  (r.m.s.)  errors 
of  the  difference  equations  over  each  computation  oyolet  to  obtain 
comparable  numbers  for  the  r.m.s.  errors,  they  have  been  normalised 
with  respect  to  the  r.m.s.  change  in  aircraft  pitoh  rate  from  its 
value  at  the  begining  of  the  relevant  computation  cyole. 

At  the  end  of  the  first  three  computation  oyoles  the  r.m.s. 
error  of  the  z-differenoe  equation  was  less  vhan  that  of  the  other 
two,  so  that  AKv  was  reduoed  to  4.2  dB  in  three  steps  of  0.6  dB. 
However,  at  the  end  of  the  fourth  cycle  the  y-differenoe  equation 
had  thp  smallest  error,  and  consequently  no  gain  change  was  made. 

It  may  be  expected  that  this  will  happen  oooasionally  beoause  of  the 
approximations  made  in  the  differenoe  equations.  For  the  remaining 
three  computation  oyoles  shown  in  the  figure  further  decreases  in 
gain  were  oorreotly  generated. 

Over  the  same  period  (but  not  shown  in  the  print-out)  it 
follows  from  equation  (D.6)  that  the  values  of  elevator  effective¬ 
ness  used  in  the  differenoe  equations  (M(  ,  K{  ,  Mj  )  were  all 

x  y  s 

Increased  by  3  dB,  so  that  M|  approaohed  the  true  value  of  the 

y 
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% 

aircraft  elevator  effectivene.e.  Thus,  at  the  beginlng  of  the 
simulation  run  M,y  was  6  dB  too  .mall,  but  after  0.7  ....  it 
only  3  dB  too  small.  The  closer  matching  „f  „  to  the  true  elev¬ 
ator  effectiveness  la  reflected  In  greater  accuracy  m  the  dlfi.r- 
ence  equation  solutions.  For  erample,  at  the  end  cf  the  first 

r^TnooTl1’  i"dlffer,n0*  —**«  »~r  was  49  per  cent 
(0.0488/0.0991),  whereas  at  the  end  cf  the  seventh  cycle  it  ..a  only 

0.76  per  cent  (0.0030/0.3926),  both  the  percentage  error,  with 

respect  to  the  aircraft  pitch  rate,  and  the  absolute  error  were 
reduced. 


TIME 


o.i 

AIRCRAFT  FITCH  RATE 

0.0000 

0.0019 

o.oioo 

0.0256 

0.0< 

X  DIFF.  EON.  ERROR 

Y  DIF:.  EON.  ERROR 

0. 

0. 

0. 

0. 

-0.0047 

-0.0030 

-0.0152 

-0*0103 

—0.01 

-0.0J 

Z  DIFF.  EQN.  ERROR 

0. 

0. 

-C.OOlf 

-0.0054 

-0.0) 

0«  2 

AIRCRAFT  PITCH  RATE 

0.1022 

0.1054 

0.1093 

0.1137 

0.11 

X  DIFF.  EQN.  ERROR 

0. 

0. 

-0.0004 

-0.0013 

-0.0( 

Y  DIFF.  EQN.  ERROR 

0. 

0. 

-0.0003 

-0.0008 

-o.ot 

Z  DIFF.  EON.  ERROR 

0. 

0. 

-0.0001 

-0.0003 

-0.01 

0.3 

AIRCRAFT  PITCH  RATE 

0.1663 

0.1768 

0.1875 

0.1979 

0.2C 

X  DIPF.  EQN.  ERROR 

0. 

0. 

-0.0001 

0.0001 

o.ot 

Y  DIFF.  EQN.  ERROR 

0. 

0. 

-0.0001 

-0.0000 

o.ot 

l  DIFF.  EQN.  ERROR 

0. 

0. 

-0.0001 

-0.0001 

-o.ot 

0.4 

AIRCRAFT  PITCH  RATE 

0.2403 

0.2415 

0.2423 

0.2429 

0.24 

X  DIFF.  EQN.  ERROR 

0. 

0. 

0.0004 

0.0009 

o.oc 

Y  DIFF.  EQN.  ERROR 

0. 

0. 

0.0003 

0.0007 

o.ot 

Z  DIFF.  EON.  ERROR 

0. 

0. 

0.0001 

0.0004 

o.ot 

0.5 

AIRCRAFT  PITCH  RATE 

0.2683 

0.2763 

0.2850 

0.2938 

0.3C 

X  DIFF.  EQN.  ERROR 

0. 

0. 

-0.0004 

-0.0011 

-o.oc 

V  DIFF.  EQN.  ERROR 

0. 

0. 

-0.0003 

-0.0007 

-o.oc 

Z  DIFF.  EQN.  ERROR 

0. 

0. 

-0.0001 

-0.0003 

-o.oc 

0.6 

AIRCRAFT  PITCH  RATE 

0.3323 

0.3326 

0.3321 

0.3310 

0.32 

X  DIFF.  EQN.  ERROR 

0. 

0. 

0.0006 

0.0017 

O.OC 

Y  DIFF.  EQN.  ERROR 

0. 

0. 

0.0004 

0.0010 

o.oc 

Z  DIFF.  EQN.  ERROR 

O.T 

0. 

0. 

0.0001 

0.0003 

o.oc 

AIRCRAFT  PITCH  RATE 

0.3367 

0.3421 

0.3484 

0.3555 

0.3< 

X  DIFF.  EQN.  ERROR 

0. 

0. 

-0.0006 

-0.0017 

-O.OC 

V  DIFF.  EQN.  ERROR 

0. 

0. 

-0.0003 

-0.0008 

-o.oc 

Z  DIFF.  EQN.  ERROR 

0. 

Fig. 

0. 

A.l 

0.0000 

0.0001 

o.oc 

Sample  cf  the  Simulation  Program  Print-Out 


T 

r 


RHSX 

RMSY 

RMSZ  GAIN 

ERRDR(DB) 

0. 

6824  0. 

4963  0. 

311) 

6.0 

!  0.0019 

0.0100 

0.0256 

0.0453 

0.0643 

0.0793 

0.0893 

0.0954 

0.0991 

1  °* 

-0.0047 

-0.0152 

-0.0292 

-0.0427 

-0.0534 

-0.0602 

-0.0639 

-0.0658 

0. 

-0.0030 

-0.0103 

-0.0203 

-0.0303 

-0.0364 

-0.0438 

-0.0470 

-0.0488 

0. 

-0.0014 

-0.0054 

-0.0114 

-0.0178 

-0.0234 

-0.0274 

-0.0301 

-0.0319 

0. 

3121  0. 

1900  0. 

0722 

5.4 

0.1054 

0.1093 

0.1137 

0.1187 

0.1244 

0. 1308 

0.1382 

0.1467 

0.1561 

0. 

-0.0004 

-0.0013 

-0.0025 

-0.0042 

-0.0065 

-0.0095 

-0.0132 

-0.0176 

0. 

-0.0003 

-0.0008 

-0.0015 

-0.0026 

-0.0040 

-0.0058 

-0.0080 

-0.0107 

0. 

-0.0001 

-0.0003 

-0.0006 

-0.0010 

-0.0015 

-0.0022 

-0.0030 

-0.0041 

0. 

1754  O.l 

0992  O.i 

0197 

4.6 

0.1768 

0.1875 

0.1979 

0.2076 

0.2164 

0.2239 

0.2301 

0.2348 

0.2381 

0. 

-0.0001 

0.0001 

0.0006 

0.0019 

0.0040 

0.0071 

0.0113 

0.0165 

0. 

-0.0001 

-0.0000 

0.0003 

0.0009 

0.0021 

0.0039 

0.0063 

0.0094 

0. 

-0.0001 

-0.0001 

-0,0002 

-0.0001 

0.0001 

0.0005 

0.0012 

0.0021 

o.; 

2340  0.1 

1042  0. 

1633 

4.2 

0.2415 

0.2423 

0.2429 

0.2438 

0.2451 

0.2473 

0.2507 

0.2553 

0.2612 

0. 

0.0004 

0.0009 

0.C014 

0.0015 

0.0010 

-0.0002 

-0.0023 

-0.0053 

0. 

0.0003 

0.0007 

0.0011 

,  0.0014 

0.0014 

0.0010 

0.0001 

-0.0012 

0. 

0.0001 

0.0004 

0.0007 

0.0011 

0.0015 

0.0019 

0.0022 

0.0024 

0.0692  0.0337  0.0223 

4.2 

0.2763 

0.2850 

0.2938 

0.3025 

0.3105 

0.3176 

0.3235 

0.3279 

0.3308 

0. 

-0.0004 

-0.0011 

-0.0015 

-0.0016 

-0.0010 

0.0004 

0.0029 

0.0064 

0. 

-0.0003 

-0.0007 

-0.0010 

-0.0012 

-0.0009 

-0.0002 

0.0010 

0.0029 

0. 

-0.0001 

-0.0003 

-0.0006 

-0.0008 

-0.0010 

-0.0012 

-0.0012 

-0.0010 

1.2100  0.9569  0.6846 

3.6 

0.3326 

0.3321 

0.3310 

0.3296 

0.3285 

0.3278 

0.3281 

0.3296 

0.3324 

0. 

0.0006 

0.0017 

0.0029 

0.0041 

0.0050 

0.0053 

0.0049 

0.0035 

9. 

0.0004 

0.0010 

0.0019 

0.0027 

0.0035 

0.0040 

0.0042 

0.0040 

0. 

0.0001 

0.0003 

0.0007 

0.0012 

0.0018 

0.0025 

0.0032 

0.0040 

• 

0.1310  0.0706  0.0192 

3.0 

0.3421 

0.3484 

0.3555 

0.3630 

0.3705 

0.3776 

0.3838 

0.3889 

0.3926 

0. 

-0.0006 

-0.0017 

-0.0030 

-0.0044 

-0.0055 

-0.0061 

-0.0059 

-0. 0048 

0* 

-0.0003 

-0.0008 

-0.0015 

-0.0022 

-0.0028 

-0.0032 

-0.0033 

-0.0030 

0.0000 

0.0001 

0.0000 

-0.0001 

-0.0003 

-0.0006 

-0.0010 

-0.0014 

■  Program  Print-Out 
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«•*»  ADAPTIVE  control  system  simulation  program  ***» 


#  MIMIC  SECTION 

fSETUP  15  DISK. 1496 

SDATA 

PAR I WVEL  *UAC . TS ) 

PAR  IGF. GV, STOP. K. TEMPO) 
PARIMQ.MALF .MDEL.LALF .LDEL ) 

DT  «  0.002 

DTMIN  «  TS 
DTMAX  *  TS 
WIND  GUST  GENERATOR 

TC  *  666. /U AC 

SD  ■  SORUOO.O*!  1.-EXPI-2.FITS/TC)  H  ) 

SNO  «  145.91 

RAN  «  SR2 1 0. . SD.SNO. T ) 

WVEL  *  RAN+WV£L*EXPt-TS/TC> 

ALFN  «  -WVEL/UAC 
PILOT  INPUT 

QC  ■  0.0 

MODEL 

OMOD  ■  FTRIOC. .5  I 
ERROR  SUMMATION 

EPSY  »  QMOD-QGYRO*AMPL*SINt30.0*T) 
COMPENSATOR 

COMOf  ■  10.*EPSY  -9.*FTR(EPSY..005) 

FIXED  GAIN 

KFOP  *  -GF*COMOP 
VARIABLE  GAIN 

KVOP  ■  KFOP*GV 
SERVO/ACTUATOR 

DDDEL  ■  -180.#DDEL-3. 24E4* ( DEL-KVOP ) 

DDEL  -  INT (DDDEL. 0. ) 

DEL  -  INTIDDEL.O. ) 

EL  «  INTI  DEL  *0. ) 

AIRCRAFT  SHOFT-PERIOD  DYNAMICS 

DQAC  •  MQ*QAC*MALF#ALF+MDEL#EL 

OAC  *  INTIDOAC.O.) 

DALF  -  QAC-LALF*ALF-LDEL*EL 

ALF  ■  INTIDAlF.O. )+ALFN*57.3 

DGAMA  *  OAC-DALF 

NZG  *  ( DGAMA*UAC > / ( 57.3*32.2 ) 

RATE  GYRO 

DDGYRO  -  -40O.*DGYRO-1.6E5#(QGYRO-QACI 
DGYRC  -  INT I DDGYRO. 0. ) 

OGYRO  ■  INTIDGYRO.O. ) 

E-FILTER 

E  «  EL 

O-FILTER 

Q  -  OGYRO 


53 


UA/EE/67-2 


MI MT C-rORTRAN  COMMUNICATION 

GV  SRl(E.Q.GV.MDEL.TEMPG.K) 

GVERR  a  TEMPG+T-T 

F  IN  <  T  *  STOP ) 

OUTPUT  TO  PLOTTER 

PLO< T. GVERR. EL. GAC.NZG) 

ZER<5. 0*6. 0.9.0* 12.0* 15.0) 

SCA (0.5 .5.0,10.0.0. 1.0. 1  ) 

PLO ( T  *  ALF  » WVEL ) 

ZERf  5.0*  28.0*21.0) 

SCA( 0.5,5.0.40.0) 

END 

C**»*  FORTRAN  SECTION 
$ IBFTC  MM  1 3 

SUBROUTINE  SRI <A*B*C*D»GVERR*ADJGVK*GV) 
DIMENSION  P ( 95  5 ♦ R ( 2500 ) »S ( 2500 ) #FF ( 9100 » 

COMMON  P • R • S , FF  » I  OUT  *  I  PAR  » I NOUT , I HDR  *  I F I N  •  I  END  » 
C  NPAR*LKDR*LKDR2*1LC.IPC.IS*IC0UNT*ISKIP 

DIMENSION  EUO)  *Q(  10)  *X<10)  *YU0>  *21 10) 
DIMENSION  XERR(IO) *YERR< 10) *ZERR( 10) 

DIMENSION  EA ( 50 ) *  0BI50) 

DATA  1/0/ 

IFUOUT.NE.l  )  RETURN 
IFfR(l). NE . 0. )  GO  TO  50 
I  »0 

CMDELo-D 
G  V*C 

REQDGV*52.95/CMDEL 
GVRAT*GV/REQDGV 
I F ( GVRAT  —  1 . )  1.2.3 
1  GVERR*  -8.686*ALOG< l./GVRAT) 

GO  TO  4 
?  GVERR*0. 

GO  TO  4 

3  GVERR* 8 . 686* A LOG  <  GVRAT ) 

4  K*AOJGVK 

WR I TE ( 6 , 20 ) 

20  FORMAT!/. 5H  TIME. 75X .4HRMSX .6X .4HRMSY.6X .4HRMSZ  » 
C  16M  GAIN  ERROR ( DB ) *//  ) 

RETURN 
50  I»I+1 

C**#*  RECEIVE  DATA  FROM  MIMIC  SECTION 
EA ( I  )«A 
OBf I >*B 

IFC I .NE.50 ) RETURN 

r#*##  SAMPLES  FOR  THE  DIFFERENCE  EONS. 

J*  0 

DO  70  L ■ 1 ♦ 46 . 5 
J*J+1 

E< J)-EA(L) 

70  Q  <  J )  *OB  ( L  ) 
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C****  COMPUTED  VALUES  OF  AIRCRAFT  CONSTANTS 

!SUMSDO=0« 

SUM5QX*0 • 

SUMSOY*0. 

.  SUMSQZ®0* 

DMY*S2.95/GV 

DMX=DMY/2.0 

DMZ*DMY*1.5 

Df$MY*8» 686*ALOG  (  DMY/ *  2193  ) 

D8NX*DBMY-6« 

08MZ* DBM Y+3 • 522 
IF(GV.GT.120.7)DBMX«0*0 
IF(GV.LT*1*5 )DBMZ»47.66 
ZWAX*«002#FXP ( • 1 5 18*DBMX ) 

Z WAY **00 2* EX Pi •1518+DBMY ) 

Z VAZ*,002*EXPl ,1518+DBMZ  > 

WAX«OMEGA ' DBMX ) 

WAY*OMEGA ( DBMY ) 

WAZ*0ME6A ( DBMZ ) 

C#**«  COEFFICIENTS  FOR  THE  CENTRAL  DIFFERENCE  EQUATIONS 
DMX=-DMX 
DMY*~DM Y 
DMZw-DMZ 

CX*l./( 1#+.01*ZWAX> 

CXE  31  * • 005*DMX*CX 
CXE2*.0001*DMX*ZWAX*CX 
CX2*(2.-.0001*WAX*WAX)*CX 
CX1*( 1»0-0*01*ZWAX>*CX 

c 

CY«  l*/( 1*+#01*ZWAY) 

|  CYE31».005*DMY*CY 

j  CYE2»«0001*DMY*ZWAY*CY 

1  CY2»(2*-.0001*WAY*WAY)*CY 

!  CY1-C 1*0“0.01«ZWAY)*CY 

i  C 

CZ*l./< 1.+.01+ZWAZ) 

j.  CZE3l«.005+DMZ*CZ 

CZ62*,0001*DMZ*ZWAZ*CZ 
CZ2«(2,-,0001*WAZ*WAZI»C Z 
CZ1*< 1*0-0#01*ZWAZ)#CZ 

!  Ct***  ITERATION  OF  THE  DIFFERENCE  EQUATIONS 

!  DO  100  J«1.2 

I  XU)-O(J) 

Y(J)«X(J) 

100  Z(J}*Y(J) 

DO  200  J»lt8 

XI J+2)«CX£31*(El J+2)-E( J) )+CXE2*Ef J+1>+CX2*X( J+l)-CXi*X( J) 

Y( J+2)-CYE31*(E< J+2)-E( Jt >+CVE2*E{ J+l )+CY2*Y( J+l )-CYl*YC J) 

H  J+2)*CZE31*(E<  J+2)-E<  J)  >+CZE2*E< J+l  >+CZ2*Z(  J+l  >-CZI«Zl J) 

SUM$QX«SUMSQX+ l X ( J+2 I -Q ( J+2 ) ) **2 

iUf ’SOY * SUMSQ Y +  l  Y  l  J+2 ! -0 ( J+2 )  I  **2 

SUMSQZ-SUMSQZ+ ( Z l J+2 ) -Q ( J+2 )) **2 

$UMSD0*5UMSDQ+<Q( J+2)-Q(2) )**2 
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CALCULATION  OF  RMS  ERRORS 
RMSX-SQRT ( SUMSQX/SUMSDQ ) 

RMSY-SQRT ( SUMSQY/SUMSDQ) 

RMSZ-SQRT  < SUMSOZ/SUMSDQ) 

C*#**  DECISION  W' IE  THE  R  TO  CHANGE  GAIN 
PRFRR-GVl-RU 
I F { K  —  1 ) 16*30*16 

30  IF( (RMSX.GT.RMSY).AND. (RMSY.GT.RMSZ) )GO  TO  9 
I F ( ( RMSX.LT.RMSY ) • AND* ( RMSY.LT.RMSZ ) I GO  TO  11 
GO  TO  16 

C*#**  GAIN  WAS  TOO  HIGH 

9  GV-GV/1.072 

I F ( GV-1 •  )  10*8.8 

10  GV-1. 

GO  TO  8 

11  IFIRMSZ.LT. ( 3.0*RMSX) )GO  TO  16 
C*##*  GAIN  WAS  TOO  LOW 


GV«GV*1.035 

IFIGV-2'  .1  .4)  8.8*12 

12 

GV-241.4 

8 

3VRAT-GV/REQDGV 

I F ( GVRAT-1 • )  13,14*15 

13 

GVERR*  -8.686*AL0G( l./GVRAT) 

GO  TO  16 

14 

GVERR-O. 

GO  TO  16 

15 

GVERR-8. 686*ALOG ( GVRAT ) 

c*##* 

OUTPUTS  TO  BE  PRINTED 

16 

DO  300  M-1.10 

Xi£RR(M)=X(M)  -Q  (M ) 

YERR(M)*Y(M) -0 (M ) 

300 

ZeRR(M)»Z(M)-Q(M) 

WR I TE ( 6  *22 ) R  < 1 ) ♦ RMSX .RMSY .RMSZ 

*PRERR 

22 

FORMAT (F4.1.70X.3F10. 4. F16.1) 
WRITE (6. 60) (Q(M) , M-1.10) 

60 

FORMAT ( 20H  AIRCRAFT  PITCH  RATE 
WRITE (6 *61 ) (XERR(M) *M«1»10) 

•10F10.4. 

61 

FORMAT ( 20H  X  DIFF.  EON.  ERROR 
WRITE (6. 62) (YERR(M) *M-1*10) 

* 10F10.4) 

62 

FORMAT ( 20H  Y  DIFF.  EON.  ERROR 
WRITE (6.63 ) (ZERR(M) *M-1.10) 

.10F10.4) 

63 

FORMAT ( 20H  Z  DIFF.  EON.  ERROR 
1-0 

RETURN 

END 

.10F10.4.//) 

SIBFTC  OME 

FUNCTION  OMEGA (DBGEN) 
IPIDBGEN.GT.il.)  GO  TO  53 
OMEGA- .5*,03*DBGEN 

RETURN 

53 

I F ( 45. 5.GE .DBGEN )  GO  TO  56 
OMEGA-DBGEN-40.6 

RETURN 

5* 
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56  OMEGA*-. 57+,121*DBGEN 
RETURN 
END 

SUBROUTINE  SR2  ( XMEAN»STD  * ISTART  *D»E*F*ANS) 
I  START- 1 START*305 17578125 
ANS-OR(ARS( I  START *8 ) »2«*34> 

CALL  MI MRN ( XMEAN * STD» ANS ) 

RETURN 

END 
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Appendix  B 

Digital  Simulation  of  the  Effect  cf  Vertical  Random  Wind 
Ousts  on  the  Performance  of  the  Gain  Computer 

Effect  of  Gusts  on  the  Airframe 

For  flight  through  atmospheric  turhulenoe  only  the  vertioal 

components  of  wind  velocity  (W  )  have  any  influenoe  on  the  aircraft 

8 

longitudinal  motion.  The  principle  effect  is  that  a  downdraft 
(positive  in  the  NACA  axis  convention)  results  in  a  negative  incre¬ 
mental  change  (o  )  in  the  angle  of  attaok,  owing  to  a  change  in 

o 

direction  of  the  relative  wind. 

Wg 

“g  *  “  57.3  ^  (deg.)  (B.1) 

where  U  is  the  aircraft  forward  velooity. 

Spectral  Density  of  Turbulanoe 

Over  a  short  period  of  time  variations  in  atmOBpherio  turbul¬ 
ence  may  be  assumed  to  have  a  stationary  Gaussian  distribution 
within  a  local  region  of  the  atmosphere.  This  assumption  implies 
that  the  region  or  time  is  small  relative  to  the  entire  flight  path 
or  flight  duration,  but  is  large  enough  for  statistical  equilibrium 
to  he  achieved  (Ref.  18*5).  At  any  instant  within  suoh  a  region 
along  the  flight  path  the  vertioal  components  of  wind  velooity  are 
characterized  by  the  root-mean-equare  (r.m.s.)  intensity,  and  by 
the  spatial  distribution  of  velooity.  By  Fourier  analysis,  this 
spatial  distribution  may  be  decomposed  into  sine  waves,  eaoh  of 
whioh  has  an  associated  wave  number.  The  predominant  wavelength, 
the  one  which  contains  the  largest  turbulenoe  energy,  has  typioally 
a  length  (  Ap)  of  0.8  miles,  so  that  the  frequency  (f  )  st  whioh  an 
airoraft  traverses  these  waves  is, 


An  approximate  expression  (Ref.  15»2)  for  gust  power  spsotral 
density  (  $  (w),  ft?/sso?  per  radian/seo. )  whioh  rsflseis  the  influ- 
enoe  of  flight  speed  on  the  received  distribution  of  gust  velooity 
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is  given  below,  _ 

4  T  W* 

0J")  -  - *~lTi  for  0< w  <  «o  (B.3) 

s  1  +  w*t! 

o 

where,  w  »  2nf  (rad. /see.) 

2  2  2 

Wt.  «  mean-square  gust  velooity  (ft. /sec.) 

O 

and,  Tg  »  (see.,  for  U  in  ft. /sec.) 

Equation  (B.3)  has  been  used  by  Cornell  Aeronautical  Laboratory 
(Ref.  3*25)  for  testing  the  gust  response  of  a  proposed  supersonic 
transport  aircraft,  and  is  the  expression  for  spectral  density  which 
was  used  in  this  simulation. 

Analog  Simulation  of  Gusts 

Before  the  digital  simulation  of  gusts  is  described  the  more 
familiar  problem  of  analog  simulation  will  be  mentioned,  because  it 
illustrates  the  relative  simplicity  of  all-digital  simulation. 
Equation  (B.3)  may  be  written, 

**<’>  •  4  Ts  *1  |rrj^|2  o-4> 

Prom  the  form  of  equation  (B.4)  it  is  clear  that  the  required 
spectral  distribution  could  be  obtained  by  analog  means  by  passing 
white  noise  of  (uniform)  power  density  4Tgw|  through  a  first  order 
filter  of  time  constant  Tg.  Although  the  bandwidth  of  praotical 
white  noise  generators  is  neoessarily  restricted  (e.g.  DC  to  40 
cyoles/sec,  for  Model  401  from  Electronic  Associates  Inc.),  the 
filter  output  approximates  the  required  noise  characteristics  suffic¬ 
iently  accurately  for  most  purposes  as  long  as  the  noise  generator 
has  a  uniform  power  output  up  to  a  frequency  of  a  decade  or  more 
above  the  filter  cut-off  frequsnoy. 

One  method  of  obtaining  samples  of  band-limited  noise  for  use 
in  a  digital  oomputer  program  is  to  reeord  the  filter  output  on 
magnetic  tape  via  an  analog-to-digital  converter,  but  there  are 
several  dieadvantagee  to  this  techniques 

a.  The  sampl 4ng  time  intervals  (t  )  used  in  making  the  tape 

B 
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determine  the  integration  step  interval  whioh  must  be  used  in 
the  program,  and  hence  the  aoouracy  of  the  simulation.  Any 
difference  between  the  two  time  soales  alters  the  effective 
cut-off  frequency  of  the  filter,  and  so  changes  the  noise 
characteristics. 

b.  A  new  tape  is  required  for  eaoh  Plight  Condition,  since 
the  filter  cut-off  frequency  is  a  function  of  the  airoraft 
velocity. 

c.  There  is  a  physical  problem  involved  in  making  the  tapes 
and  ensuring  compatibility  with  the  computer  installation. 

Digital  Simulation  of  Gusts 

A  simple  method  of  generating  band-limited  noise  samples  by 
means  of  program  statements  is  to  employ  the  Gaussian  random  number 
subroutine.  The  method  outlined  below  ensures  that  the  resulting 
wind  gust  samples  (Wg^,  n  *  0,  1,  2,  ...)  have  the  required  spectral 
density,  as  specified  by  equation  ( B. 3). 

Although  the  physically  measurable  wind  gust  spectral  density 
exists  only  over  the  range  of  positive  frequencies,  it  is  usual  and 
convenient  for  the  purpose  of  mathematical  analysis  to  define  a 
spectral  density  function  0'(w),  such  that. 


0'(w) 


V2  *g(") 


for  04  w<  «o 


l/2  0_(-w)  for  -ee<w^0 


(B.5) 


so  that  (w)  is  an  even  funotion. 

The  autocorrelation  function  for  0' (w)  is  defined  as. 


VT>  ■ 


1  t.  tf  f atisa 

**  Ml 


(B.6) 


Henoe,  by  the  use  of  oontour  integration. 


Ur) 


if  exp(-hVTg) 


(B.7) 


1 


G A/EE/67- 2 


I 

I 


In  equation  (B. 7)  the  absolute  value  of  X  has  been  taken 
beoause  the  autocorrelation  function  is  always  even*  It  may  be 
noted  that  the  value  of  the  autocorrelation  function  for  X  «  0 
yields  the  mean-square  gust  velocity.  A  further  useful  property 
of  the  autocorrelation  function,  which  will  now  be  used,  derives 
from  the  faot  that  all  waveforms  which  have  the  same  power  spectral 
density  distribution  also  have  the  same  autocorrelation  funotion 
(Ref.  11*58). 


Suppose  that  the  wind  gust  samples  are  formed  at  tiiae  intervals 
t0  from  the  sequenoe  of  random  numbers  *n  (n  »  0,  1,  2,  ..*),  having 
a  specified  varience  a 2  and  sero  mean,  in  the  following  way* 


®n+l 


*„♦!  ♦  VIP(‘P*«) 


(B.8) 


where  the  term  exp(-0t8)  represents  a  decay  factor  for  the 

previous  gust  sample,  and  the  starting  value  I-  may  be  arbitrarily 

°o 

chosen.  It  may  be  shown  (Ref.  7*1-130, 1-131 )  that  at  the  disorete 
delay  times  given  by  T  «  nt  the  autocorrelation  funotion  for  the 

13 


samples  Wg^  Is, 


R(r) 


1  -  exp(-2pt8) 


(B.9) 


Henoe,  if  R(?)  is  to  be  Identical  to  Rg(T)  it  is  neoessary  that, 


exp(-irf/Tg) 
Thus,  P 


x2  gjn?(-fllT( ) 
1  -  exp(-2pt0) 

lAg 


(B.10) 


and, 


"  »|fl  -  exp(-2t8/Tg)]  (B.ll) 


Use  of  the  find  Oust  Equations 

To  test  the  performance  of  the  gain  oomputer  in  the  presenoe  of 
random  inputs  a  wind  gust  intensity  of  r.m.s.  value  20  ft./seo.  was 
seleoted,  einoe  this  magnitude  was  used  by  Honeywell  (Ref.  9*4)  in 
the  development  of  adaptive  oontrol  systems  for  the  X-15  end  X-20 
(Dyne  Soar)  vehloles.  For  the  integration  step  Interval  (t8)  used 
in  the  MIMIC  simulation  program  a  slse  of  0.0005  eeo.  was  ohosen  as 
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a  compromise  between  accuracy  of  computation  and  time  taken  by  the 
IBM  7094  computer  to  complete  a  simulation  run.  Then,  for  each 
Plight  Condition  the  time  constant  T^  was  calculated  from  the 
aircraft  velocity,  as  shown  under  equation  (B. 3).  Substitution  of 
the  values  for  t  ,  and  T  into  equation  (B.ll)  gave  the  varience 
required  for  the  random  number  subroutine.  The  wind  gust  samples 
were  then  calculated  from  the  random  numbers  by  the  use  of  equation 
(B.8).  Finally,  samples  of  the  disturbance  angle  of  attack  were 
found  by  substituting  the  wind  gust  values  into  equation  (B.l). 

Results  of  the  Simulation 

Figures  B.l  through  B. 4  show  the  aircraft  response  and  perfor¬ 
mance  of  the  gain  computer  in  the  presence  of  random  wind  gusts. 

The  dither  signal  was  used,  but  there  were  no  pilot  command  inputs. 

It  should  be  noted  that  the  vertioal  soales  in  Fig.  B.4  are  differ¬ 
ent  from  those  in  the  other  three  figures,  and  in  all  four  the  soale 
of  4Ky  ic  twioe  that  which  was  earlier  used,  so  that  the  finite 
step  ohanges  in  gain  are  visible.  Starting  from  the  top  trace,  the 
curves  are* 

a.  Wind  gust  velocity,  Wg. 

b.  Incremental  change  in  angle  of  attack,  a. 

o.  Normal  acceleration,  n  . 

z 

d.  Aircraft  pitoh  rate,  • 

e.  Elevator  defleotion  angle,  $  . 

f.  Control  loop  gain  error,  AK^. 

Aircraft  gust  sensitivity,  in  terme  of  its  normal  aooeleration 
and  pitoh  rate  responses,  is  known  (Rsf.  14:13)  to  be  primarily 
dependent  on  the  lift  ourve  elope  (proportional  to  LQ)  and  the 
etatio  stability  (proportional  to  Ma),  reepeotively.  These  relation¬ 
ships  srs  apparent  when  the  and  ourves  Are  compared  with  the 
magnitudes  of  the  Lq  and  stability  derivatives  (Appendix  C, 

Table  I)  for  eaoh  of  the  Flight  Conditions.  For  example,  in  Flight 
Condition  28  (Fig.  B.l),  where  L  takes  on  its  largest  value,  the 
normal  aooeleration  response  is  also  the  greatest  of  the  four,  with 
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peaks  of  ±2.5g.  For  a  manned  airoraft  this  acceleration  level  is 
obviously  unacceptable,  and  in  practise  is  unlikely  to  arise, 
because! 

a.  An  r.m.s.  gust  intensity  of  20  ft. /sec.  represents  severe 
atmospheric  turbulence,  for  vhioh  the  probability  that  it  will 
be  encountered  is  less  than  0.01  per  cent  (fief.  15*3). 

b.  In  simulating  the  airoraft  dynamics  it  wets  assumed  that 
the  &  stability  derivatives  were  negligibly  small  (Appendix  C). 
This  assumption  implies  that  the  pressure  distribution  over  the 
wing  changes  instantaneously  when  the  angle  of  attaok  is 
suddenly  changed  (Ref.  5*158).  For  pilot  commands  transmitted 
through  the  pre-filter  model  reference  the  rate  of  change  of 
angle  of  attack  is  sufficiently  slow  for  the  assumption  to  be 
valid.  However,  as  the  Figures  show,  wind  gusts  can  cause  very 
rapid  changes  in  the  angle  of  attaok,  and  the  assumption  is 
then  no  longer  strictly  correot.  Qualitatively,  if  the  aero¬ 
dynamic  lag  is  taken  into  consideration,  it  may  be  seen  that 
the  peaks  on  the  curves  in  Figs.  B.l  through  B*4  would  be 
considerably  smoothed. 

To  ensure  that  the  airoraft  and  control  system  had  achieved  a 
'steady-state*  all  the  simulation  runs  with  wind  gust  inputs  were 
made  15  eeo.  long,  of  whioh  the  Figures  show  only  the  first  6  sec. 
The  airoraft  responses  for  the  first  6  sec.  are  representative  of 
the  complete  run,  but  exoept  for  Flight  Condition  28  the  loop  gain 
decreased  further.  For  Flight  Conditions  7  end  24  the  average  gain 
error  was  approximately  -6  dB,  and  for  Flight  Condition  32  was  -4  dB. 
The  gain  for  Flight  Condition  28  oould  not  deorease  below  its  design 
value,  as  defined  by  equation  (D. 6),  since  the  variable  gain  setting 
(Ky)  was  then  on  its  lower  end  stop  for  this  oondition  of  maximum 
elevator  effectiveness.  It  may  be  noted  that  in  Figs.  B.1  through 
B.4  the  dither  signal  is  not  apparent  on  the  pitoh  rate  ourves 
beoause  of  its  small  amplitude  (Ref.  page  25). 
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Flight  Condition  28  -  Aircraft  Response  and  Performance  of 
the  Cain  Computer  for  Random  Vind  Gust  Disturbances 
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Pig.  B.2 

Flight  Condition  7  -  Airoraft  Response  and  Performance  cf 
the  Cain  Computer  for  Random  Wind  Gust  Disturbances 
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Fig.  B.3 

Plight  Condition  24  -  Aircraft  Response  and  Performance  of 
the  Cain  Computer  for  Random  Wind  Gust  Disturbances 
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Fig.  B.4 

Flight  Condition  32  -  Aircraft  Response  and  Performance  of 
the  Gain  Computer  for  Random  Wind  Oust  Disturbances 


Functions  in  Terms  of  the  Pitch  Stability  Derivatives 
The  Basic  Equations  of  Motion 

A  complete  development  (e.g.  Refs.  2  and  5)  of  the  general, 
linearised  aircraft  equations  of  motion  involves  the  use  of  a 
considerable  number  of  approximations  and  qualifying  statements, 
and  is  quite  long.  Consequently,  only  a  brief  outline  is  given, 
together  with  an  indication  of  the  principal  assumptions. 

From  the  point  of  view  of  the  control  system  designer  it  is 
convenient  if  the  equations  are  ultimately  expressed  in  transfer 
function  form.  A  complication  which  is  often  encountered  is  that 
the  equations  are  written  in  non-dimensional  form.  However,  for 
supersonic  flight  the  variations  in  magnitude  of  the  non-dimensional 
stability  derivatives  (coefficients  of  the  differential  equations) 
as  a  result  of  Mach  Number  effects  cannot  be  ignored}  in  this  event 
the  most  straightforward  approach  is  to  assooiate  with  a  number  of 
reference  Flight  Conditions  a  distinct  set  of  numerical  values  for 
the  stability  derivatives. 

The  basic  equations  of  motion  may  be  derived  in  the  following 

way* 

a.  A  body  axis  system  is  assumed  to  be  rigidly  fixed  to  the 
airframe}  the  axes  are  known  as  stability  axes  if,  when  the 
aircraft  is  flying  straight  and  level,  the  x-axis  is  direoted 
along  the  airoraft  velocity  vector,  the  y-axis  along  the 
starboard  wing,  and  the  a-axis  downwards.  The  use  of  these 
axes  enables  the  equations  to  be  derived  in  their  simplest  form 

b.  Euler’s  equations  of  motion  for  a  rotating  rigid  body  and 
Newton's  equations  for  motion  of  the  center  of  mass  are  written 
as  the  moment  and  foroe  equations,  respectively.  The  external 
foroes  aoting  on  the  airframe  are  due  to  aerodynamio  foroes  and 
gravity  foroes,  while  the  external  moments  about  the  oenter  of 
mass  are  due  to  aerodynamic  forces  alone. 
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o.  Each  of  the  foroes  and  momenta  is  assumed  to  consist  of  two 
components ;  one  a  constant  component,  determined  by  the  partic¬ 
ular  reference  Plight  Condition,  and  the  other  a  perturbation 
component  rhich  arises  when  the  aircraft  is  disturbed  from  the 
reference.  The  perturbation  values  of  the  aerodynamic  forces 
and  moments  are  assumed  to  be  linear  functions  of  the  distur¬ 
bances,  a  concept  which  gives  rise  to  the  definition  of  the 
stability  derivatives. 


d.  The  steady-state  values  of  the  forces  and  moments  are  sub¬ 
tracted  from  both  sides  of  the  resulting  equations,  and  prod¬ 
ucts  of  small  angles  and  angular  rates  are  neglected.  When 
this  is  done,  it  is  found  that  the  linearised  longitudinal 
equations  are  decoupled  from  the  lateral  equations.  A  further 


simplification  of  the  equations  results  if  the  airflow  around 
the  wing  and  tail  is  assumed  to  be  quasi-steady,  (e.g.  Pef.2*21, 


eqns.  (1-59)*  For  the  X-15  the  non-dimensional  derivative  Cz 


is  zero,  and  the  last  assumption  implies  that 
negligibly  small.) 


& 


and  C 

m. 

o 


are 


q 


e.  When  only  the  short-period  mode  is  considered  (i.e.  if  the 
comparatively  low  frequency  phugoid  mode  is  ignored)  the  small 
variations  which  occur  in  the  aircraft  forward  velocity  are 
insignificant.  Finally,  if  it  is  assumed  that  incremental 
changes  in  gravity  forces  which  result  from  small  changes  in 
pitch  attitude  are  negligible  compared  with  the  incremental 
aerodynamic  forces,  then  the  following  basic  equations  of 


motion  apply: 


(c.1) 


a  »  ^  —  L  8  —  L»S 
a  a  I  e 


(C.  2) 


where,  $  ■  pitch  rate  (deg. /sec.) 

ft 

a  ■  incremental  change  in  angle  of  attack  from 
the  reference  value  (deg.  ) 

Se  •  elevator  deflection  angle,  measured  from  the 
trim  position  (de-T. ) 


69 


/ 


GA/FE/67-2 

Values  of  the  moment  (M)  and  lift  (L)  stability  derivatives  in 
equations  (0.1 )  and  (0.2)  for  the  four  Flight  Conditions  (F.C.) 
which  were  &  idied  are  F'  in  Table  I.  Details  of  the  reference 
conditions  f.r  which  the  tability  derivatives  were  obtained  are 
listed  in  Table  II  (data  from  Kef.  12). 


Table  I 

Stability  Derivatives 


Derivative: 

Dimensions : 

% 

(sec.  ) 

II 

a 

,  -2  N 
(sec.  ) 

MS 

(sec. ^) 

L« 

(sec.1) 

(seer"1) 

F.C.  28 

-2.456 

-49*93 

■52.95 

2.527 

0.5801 

F.C.  7 

-0.2299 

-11.18 

-9.097 

0.2529 

0.04364 

F.C.  24 

-0.0342 

-3.520 

-1.741 

0.0546 

0.007145 

F.C.  32 

-0.05202 

-0.2587 

-0.2193 

0.04999 

0.01218 

Table 

II 

Reference  Flight 

;  Conditions 

Condition: 

Altitude 

Mach  No. 

Angle  of 

Velocity 

Dynamic 

Attack 

U 

Pressure 

Dimensions : 

(ft.) 

- 

(deg. ) 

(ft./eeo. ) 

(lb. /ft?) 

F.C.  28 

10,000 

1.2 

0.5 

1078 

1467 

F.C.  7 

50,000 

1.5 

5.0 

1453 

381 

F.C.  24 

100,000 

3.0 

11.4 

3014 

142 

F.C.  32 

0 

0.2 

17.0 

223 

59 

The  Short-Period  Transfer  Functions 

When  the  Laplace  transforms  of  equations  (C.l)  and  (C.2)  are 
taken  and  the  variable  a  is  eliminated,  the  transfer  funotion  whioh 
relates  pitch  rate  to  elevator  deflection  angle  is  obtainsdi 


$ 


•f 

A 
,  Vi 

| 

§ 


1 


1 

* 


* 


f 


•i 


q.(«) 

*.(■) 


1 


’  ♦ (L.  ~  *,  w 
•2*  <l«-  V*  ♦  (•".-■q1*  > 


(0.3) 


It  is  oonvenlsnt  to  write  this  equation  in  an  alternative  fora, 
in  whioh  the  airoraft  short-period  time  oonstant  (Ta)»  dsaping 


i 

I 

* 
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factor  (  X  )t  and  natural  frequency  (w  )  appear  explicitly* 

6a(s)  B  +  1/Ta 

<7^  ‘  "*  .2*2J»s*.2 

va  a  a 

where,  w^  «  -II  -  M  L  (rad?/sec?) 

7  a  a  q  a  '  7 

24fa*a  -  La"Mq  (rad. /sec.) 

and,  l/Ta  -  La«MaL</M|  (seoT1) 


(C.4) 

\ 


The  following  table  is  then  obtained  by  substituting  numerical 
values  for  the  stability  derivatives  into  the  above  equations. 


Table  III 

Fitch  Rate  Characteristics 


Parameter* 

Aerodynamic 

Natural 

Damping 

Time 

Osin  (M|) 

Frequenoy 

Factor 

Constant 

Dimensions* 

CM 

1  • 

O 

• 

a 

(rad./seo. ) 

- 

(seo. ) 

- 

F.C.  28 

-52.95 

7.492 

0.3325 

0.0720 

o. 4829 

F.C.  7 

-9.097 

3.353 

5.017 

• 

F.C.  24 

-1.741 

1.877 

0.0237 

24.91 

F.C.  32 

-0.2193 

0, 5111 

0.0997 

28.06 

A  variable  ¥  is  defined  by, 

i  -  *a  -  *  (c.5) 


The  normal  acceleration  (n^,  positive  for  eyeballs  down)  of 
the  oenter  of  mass  of  the  airoraft  is  related  to  }  by  (Ref.  20*141), 


n 


S' 


juL 

57.  ie 


(S’s) 


(C.6) 


From  equation  (C.6),  it  may  be  seen  that  the  normal  acceleration 
is  not  simply  proportional  to  pitoh  rate  while  the  angle  of  attack 
is  changing.  This  important  faot  is  emphasised  by  solving  equations 
(C.l),  (C.2),  and  (C.3)  for  the  transfer  funotion  whioh  relates 


normal  aootleration  to  pitoh  rate. 
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The  behavior  of  the  normal  acceleration  response  with  respect 
to  the  pitch  rate  response  ia  illustrated  by  considering  two  approx¬ 
imations  to  equation  (C.7)» 


a.  For  a  step  oommand  input,  the  initial  pitch  rate  response 
may  be  approximated  by  a  ramp  function.  This  assumption 
implies  that  ©a(t*)  is  a  constant,  and  henoe,  by  applying  the 
initial  value  theorem  to  equation  (0.7),  the  normal  acceler¬ 
ation  in  the  region  of  time  zero  is  obtained, 


(C.B) 


Therefore,  since  Lj  is  positive  and  Mj  is  negative,  the  initial 
sense  of  nz  is  anti-phase  to  that  of  the  change  in  pitch  rate. 
This  non-minimum  phase  phenomenon  has  a  simple  physical  explan¬ 
ation.  For  example,  to  initiate  a  climb  a  positive  pitch  rate 
is  required,  which  is  produced  by  a  negative  elevator  deflect¬ 
ion  angle.  The  immediate  result  is  a  reduction  in  the  lift 
component  generated  by  the  tail,  and  until  the  angle  of  attack 
is  increased  sufficiently  for  the  wing  to  provide  correspond¬ 
ingly  more  lift,  the  direction  of  the  normal  acceleration  is 
opposite  to  that  required. 


b. 


If  T  is  large,  there  is  a  period  of  time  after  an  initial 

ft 


transient  when  the  numerator  dynamios  associated  with  equation 
(C.7)  may  be  ignored.  In  this  event  the  equation  is  approx¬ 
imated  by, 


4«(*> 


K 


S  * 


1/T 


(C.9) 


where  K  is  a  (positive)  faotor  of  proportionality.  Thus, 
the  normal  acceleration  follows  the  pitoh  rate,  but  with  a  lagf 
only  in  the  steady-state  is  the  normal  acceleration  directly 
proportional  to  the  pitch  rate,  and  even  then  the  proportion¬ 
ality  faotor  is  a  funotion  of  the  aircraft  velooity.  It  is 
generally  acoepted  that  pilot  opinion  of  aircraft  handling 
qualities  is  strongly  influenced  by  the  form  of  the  normal 
aooeleratlon  response.  Consequently,  an  airoraft  whioh  has 
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a  pitch  rate  command  control  system  whioh  does  not  overshoot 

for  a  step  input  is  likely  to  be  rated  'sluggish'  (Ref.  20*20). 

Aircraft  Response  Curves  for  a  Step  Input 

Figures  C.l,  C.2,  and  C.3  show  the  airoraft  response  curves  for 
a  step  command  input  of  0.5  deg./seo.  For  eaoh  Flight  Condition  the 
loop  gain  was  set  at  its  design  tfalue  (defined  in  Appendix  D),  so 
that  the  phase  and  gain  margins  were  the  same  for  all  three  condit¬ 
ions.  Starting  from  the  top  trace,  the  curves  are* 

a.  Output  of  the  pre-filter  model, 

b.  Aircraft  pitch  rate,  Q  . 

ft 

o.  Difference  between  pitoh  rate  and  rate  of  change  of  angle 

of  attack,  t . 

d.  Elevator  deflection  angle,  &  . 

The  pitoh  rate  curves  show  that,  if  the  gain  computer  is  able 
to  identify  the  elevator  effectiveness  and  set  the  loop  gain 
aooordingly,  the  airoraft  will  follow  the  model  reference  with 
negligible  error.  When  the  normal  acceleration  curves  (proportional 
to  i  )  are  compared,  the  effeot  of  the  large  time  constants  (T  )  in 

ft 

Plight  Conditions  7  and  32  on  the  lag  of  acceleration  with  respeot 
to  pitch  rats  is  seen.  Also,  in  Fig.  C.3  the  non-minimum  phase 
characteristic  of  the  normal  acoeleration  response  is  apparent. 

An  inspection  of  the  numerator  of  equation  (C.4)  indicates  that 
the  pitch  rate  response  may  be  considered  to  consist  of  two  compon¬ 
ents*  one  proportional  to  the  derivative  of  the  elevator  deflection 
angle,  and  the  other  proportional  to  l/T  times  the  magnitude  of  the 

ft 

angle.  Thus,  for  values  of  signifioantly  greater  than  unity  the 
derivative  component  is  dominant,  and  this  effeot  is  particularly 
noticeable  in  Flight  Condition  32  (T,  »  28.06),  where  the  elevator 
rate  remained  almost  constant  for  about  twenty  ssoonds. 
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Appendix  D 

Design  of  a  Compensator  for  the  Control  System 
Design  Requirements 

A  block  diagram  of  the  complete  control  system  is  shown  in 
Pig.  1.,  in  which  the  pre-filter  model  reference  is  a  simple  first 
order  lag,  with  a  break  frequency  of  2  rad. /sec.  The  object  of  the 
design  was  to  make  the  overall  response,  from  pilot  command  input 
to  aircraft  pitch  rate,  follow  the  output  of  the  model  reference. 
This  may  be  achieved  by  making  the  bandwidth  of  the  closed  loop 
portion  of  the  control  system  five  to  eight  times  that  of  the  model 
for  all  Plight  Conditions. 

Control  Loop  Components 

Aircraft.  Phase  angle  -  frequency  curves  for  the  aircraft 
elevator  deflection  to  pitoh  rate  transfer  function  are  shown  in 
Pig.  D.l.  The  steep  slope  of  these  ourves  in  the  region  of  the 
aircraft  natural  frequency  is  indicative  of  the  low  damping  of  the 
airframe  without  pitch  stability  augmentation.  For  frequencies 
above  20  rad./eeo.  the  aircraft  transfer  function,  given  by  equation 
(C.4),  may  be  approximated  by, 

Mj 

*  3?  <D-l> 

Thus,  the  high  frequency  gain  of  the  aircraft  transfer  funotion  is 
directly  related  to  the  value  of  M|  •  As  a  result,  the  gain  cross¬ 
over  frequency  (where  the  loop  phase  shift  is  -180  deg.)  and  the 
gain  margin  may  be  held  constant  for  all  Plight  Conditions  by 
adjusting  the  variable  gain  element  (Ky)  inversely  to  Mj .  Further, 
since  equation  (D.l)  shows  that  the  aircraft  high  frequency  phone 
angle  is  the  same  for  all  Plight  Conditions,  the  phase  margin  will 
also  be  constant.  Hence,  the  gain  computer  ie  required  to  'measure* 
the  value  of  Mj  as  accurately  ns  possible. 

gervo/Aqtu a tor  Combination.  The  servo-cylinder  and  power 
actuator  transfer  funotions  (Hef.  12. are  shown  in  Fig.  u.2, 
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where  phase  angle  -  frequency  curves  are  plotted  for  three  values  of 
servo  gain  and  feedback.  In  each  case  the  servo  gain  was  selected 
to  give  a  servo  natural  frequency  of  180  rad. /sec.,  and  except  for 
the  case  of  unity  feedback  which  was  not  used,  the  feedback  transfer 
function  was  chosen  to  cancel  the  actuator  lag.  For  the  non¬ 
integrating  (K.j  and  Hg  in  Fig.  D.2)  and  integrating  (K^,  H^)  servos 
the  transfer  functions  for  the  servo/actuator  combination  are, 


Non- Integrating 
Servo: 


Integrating 

Servo* 


Vs> 

Vs> 

Vs) 


-  6.67 

.  L 61 


_ 31200 _ 

s 3  +  l80s2  +  32400s  +  31200 

1802 

s2  +  180s  +  1802 


(1).2) 

(».3) 


The  first  tests  on  the  closed  loop  control  system  were  carried  out 
with  the  non-integrating  servo,  but  because  a  droop  ocourred  in  the 
pitch  rate  response  (Fig.  D.7)  all  subsequent  simulation  was  done 
with  the  integrating  servo.  It  should  be  noted  that  for  the  purpose 
of  simulation  the  factor  6.67  associated  with  the  servo/ actuator 
combination  has  been  incorporated  in  the  fixed  gain  element  (Kf  ). 


Rate  Gyro.  From  Figs.  D.l  and  D.2  it  may  be  seen  that  for 
frequencies  above  20  rad./seo.  the  aircraft  and  servo/aotuator  angles 
contribute  dose  to  -180  deg. ,  or  more.  A  high  frequency  rate  gyro 
was  therefore  chosen,  so  that  it  would  have  only  a  small  phase  shift 
in  the  region  of  20  rad. /sec.  up  to  the  gain  crossover  frequency. 

The  rate  gyro  phase  angle  -  frequency  ourve  is  shown  in  Fig.  D.2, 
and  its  transfer  funotion  is  given  below, 


SnkMMWM 

*<•> 


400 2 

s2  +  400s  +  400 2 


OM) 


The  Complete  Control  Loop 


Phase  angle  -  frequency  curves  for  the  uncompensated  oontrol 
loop  (with  non-integrating  servo)  are  given  in  Fig.  D. 3»  and  in 
Fig.  D.4  the  oorrsaponding  gain  -  frequency  curves  are  presented. 
These  Figures  also  show  ths  phase  and  gain  oharaoteristioe  of  the 
chosen  compensator,  whose  transfer  funotion  is, 
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Vo(s>  _  10  (s  +  20) 

vi(s )’  *  a  4-  200 


(D.5) 


The  requirements  which  dictated  the  use  of  this  compensator 


were* 


a.  An  inspection  of  the  open  loop  phase  angle  curves  in 
Fig.  D. 3  showed  that  a  lead  network  waa  necessary  to  obtain 
the  required  bandwidth. 

b.  The  upwards  break-frequency  and  the  attenuation  factor 
had  to  be  as  small  as  possible.  This  was  to  ensure  that  high 
frequency  noise  would  not  become  troublesome. 

o.  A  lower  limit  for  the  upwards  break-frequency  was  deter¬ 
mined  by  the  resonance  peak  at  7  rad. /sec.  in  the  gain  - 
frequency  ourve  for  Flight  Condition  28.  In  the  uncompensated 

case  (Fig.  D.4)  there  was  a  15  dB  decrease  in  gain  in  the 

il ; 

ootave  between  8  and  16  rad. /sec.  By  placing  the  compensator 
upwards  breafcwfrequency  at  20  rad. /sea.  it  was  possible  to 
realise  a  13  dB  gain  roll-off  in  this  ootave  (Fig.  D.6), 
which  enabled  a  satisfactory  gain  margin  to  be  aohieved  in 
oon junction  with  a  high  loop  gain  at  the  lower  frequencies. 

The  design  value  of  the  control  loop  gain  setting  was  then 
established  in  the  following  way* 


a.  The  gain  -  frequency  curve  for  Flight  Condition  28  was 
taken  as  a  referenoe  when  both  the  variable  gain  (Kv)  and  the 
fixed  gain  (K^)  were  set  to  unity. 


b.  It  was  asstiraed  that  the  gain  computer  would  normally  be 
oapable  of  correctly  identifying  the  true  value  of  elevator 
effectiveness  for  any  Flight  Condition,  so  that  it  would  be 
able  to  set  the  gain  oorrectly  at  its  design  value. 


o*  The  variable  gain  was  changed  aa  a  function  of  the  computed 
e2«vctor  effectiveness  (M.  )  by  the  relation, 

Hj  for  F.C.  28  (-  52.95) 

Kv  r  Computed  elevator  effeoi'iveness 
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d.  Thus,  the  design  value  of  K  for  Plight  Condition  23  was 
defined  to  be  unity.  Also,  for  example,  the  design  value  for 
Plight  Condition  24  is  given  by, 

Ky  (P.C.  24)  -  52.95/1.741  -  30.42  -  29.7  dB 
where  the  true  elevator  effectiveness  for  Plight  Condition 
24  (from  Appendix  C,  Table  III)  has  been  substituted  into  the 
denominator  of  equation  (D.6).  It  follows  that,  if  the  comp¬ 
uted  elevator  effeotiveness  is  in  error  from  the  true  value, 
then  the  setting  of  Ky  will  be  correspondingly  in  error  from 
its  design  value.  Hence,  in  order  to  assess  the  accuracy  of 
the  gain  computer,  a  loop  gain  error  ( 4  Kv)  is  defined; 

K  ( actual ) 

AKv  -  20  K— d»Tg^)  (dB)  <D-7) 

e.  In  Fig.  D.4  it  can  be  seen  that  the  effect  of  increasing 
the  loop  gain  for  Plight  Condition  24  by  the  design  value  of 
Kv  (29.7  dB)  is  to  cause  the  gain  curve  for  frequencies  above 
20  rad. /sec.  to  become  00- incident  with  that  for  Flight 
Condition  28.  The  same  result  applies  to  all  other  plight 
Conditions. 

f.  The  fixed  gain  factor  was  then  set  at  a  value  of  10 
(Kf  *  20  dB),  which  provided  the  maximum  attainable  phase 
margin,  and  a  gain  margin  of  11.5  dB.  Phase  and  gain  ourves 
which  illustrate  the  net  effect  of  lead  compensation  and  gain 
adjustment  are  contained  In  Pigs.  D.5  and  D.6. 

Control  System  Bandwidth 

It  was  mention'd  jarlier  that  the  (major)  portion  of  the 
control  system  which  follows  the  model  reference  should  have  a 
dosed  loop  bandwidth  of  five  to  eight  times  that  of  the  model 
(i.e.  10  to  16  rad./seo.).  To  oonfirm  that  this  oondition  was 
satisfied,  a  closed  loop  frequency  response  calculation  was  made 
for  Flight  Condition  28,  with  the  non-integrating  servo  included 
in  the  loop.  It  was  found  that  the  closed  loop  gain  from  model 
response  (  6  )  to  airoraft  response  (  8  ),  measured  with  respeot 
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to  the  zero  frequency  gain,  was  3  dB  down  at  12  rad./eeo.  and  6  dB 
'  down  at  22  rad./eeo.,  and  was  therefore  considered  satisfactory. 

Aircraft  Response  Curves  for  a  Step  input 

When  the  non-integrating  servo  was  used,  the  oontrol  loop  was 
j  '  Type  Zer°»  Md  consequently,  a  steady-state  error  in  the  final  value 

of  the  pitoh  rate  response  was  expected.  This  error  appears  as  a 

I  dr°°P  in  the  pitoh  rat®  response,  as  shown  in  Pig.  D.7.  Since  it 

|  Wa®  "0t  Pra°tioahle  to  inorease  the  loop  gain  any  further,  an 

!  integrating  servo  was  adopted.  The  improved  response  with  the 

j  integrating  servo  is  shown  in  Pig.  D.8.  The  significance  of  the 

plots  of  the  term  i  in  these  two  Figures  is  explained  in  Appendix  C, 

|  Wh8r®  curvea  Blowing  the  step  responses  for  the  other  three  Plight 

j  Conditions  are  contained. 

J 
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s  lo  It  a*  #•  *0  («  loo  Tad.. 


S«C. 

Pi*.  D.5 

Phase  Angle  -  Frequenoy  Curves  for  the  Compensated 

Open  Loop 


Fig.  D.6 

Gain  -  Frequenoy  Curves  for  the  Conpeneated 
Open  Loop  With  at  its  Design  Value 


86 


GA/EE/67-2 


1  — 

WMmi 

Iaaaa* *##* 
iltlHflll 

:s»! 

!fflSnR! 

aaaaa aaaa 

nnf 

::::: 

ill 

■si iii: 

Ijipnp 

i^Ra^^R 

!HSS::s?s 

!^!5T 

lira 

::::::::: 

i«3**liMa 

IS 

P 

i iiiii ill iiiii IS 

ifl/iSE:  Hi 

!*•*•• ***■ 

•  MMIIMI 
■  •MMMM 

aaaa 
aaaa 
•  •aa 

::::::::: 

ill 

iii 

■  aa 

SSHSI 

at*  aaaa 

;*>•■§« 

mm 

:ui3:u: 

1 aaaaa aaa* 
•■*•** aiaa 

i*i***iia**i 
■  •  i*aa*  ***** 

iikaaiiaik 

:355j 

aaaSsu 

!|!l! 

:i:5j!i!:ii:::i! 

9* 

BO 

HHK: ih: 
:::: 

■aaaia Kaaa 
taliMMII 

nialitBia 

»  a ■»* •*•* 

’■•■•■■■St 

ItIHI  #BB* 

ijjijjjHj 

ill* 

HI'. 

i:k 

■hj 

cttr 

tlssilt:*: 

::::::::: 

HHl!::: 

::::: 

Saata 

;::j: 

iibs 

■\l:V 

IU 

m 

li 

i 

!«■ task 

555  Hi: 

::::::: 

:::ks: 

aa*  aaaa 
aa* aaaa 
■  Ea  kaaa 

■aaaa*** aa 
kaaaa*  aaaa 

'■aatfg  aaaa 
aaaaiaaaa 
aaaaa  kb aa 
aaaaa  aaaa 
aaaa*  aaaa 
aaakaaaa* 
aaaa*  «aaa 

shnssts 

!sjn:::» 

:::::::::: 

u» 

IHH1SP 

lipiijl 

jpiHi 

m 

iaiaaaaaaaa* 

illililiil; 

mm 

>■***■  ■•*■ 

SSaasaiaai 

BsinrsMi 

*a*«*a*al 

t:in  sms: 

a  kaiai 

iii 

Iii 

:  Hj 

iaaaa 

ra 

iiiii 

ill 

S:s::s:r:::::::: 

•**an*I**iik**ii 

aaaaa* i»**tik**l 
iaaa**-iaaii*iaai 

aaaa 

Saaa 

ssfi 

saaaaatat 
aaaaa aaaa 

tints::: 

:::: 

aaaa 

ES»i:g 

aaih  *8tt 

ill 

talSIkll**! 

5H:l!H5li 

3555  !53i 

{aSaai 

:5i::i 

in:: 

<  aaaa 
i  aaaa 
iaaaa 

:i|: 

n^S.r::3: 

Saaa 

aaaa 

jaaa 

tnt::::t 

aata 

::::: 

n::z::n 

aaaa*  a* a* 
::::::::: 

::::::::: 
aaaaa aaaa 
■aaaa aaaa 

::!5?H5: 

ill 

iill! 

H|5:ii!5i 

kaaaa* 

tasaai 

IH55! 

!E:| 

isnt 

iaaBa*  aa** •••*■• 
■aSaaa ■■•kiitiki 

I:  ' 

:::::H35: 

aaaa 

355 

aaaaa  aaaa 

aaaaiaaaa 
Kasai aaaa 
•a jaalaaa 

a  aaaa 

11:11 

iEEjjir 

iHfiiiiii 

111 

PIP 

111! 

ill 

itsss: 

EH 

k aaaaa a t* ■*£■•** 
■*■■■■ ■■*■**** *■ 
:::::::::: 

J;j| 

aaaa' 

Salaaaaaa 
aaaaa  aaaa 

!i!H 

jiiliiilj 

itsitlflf 

nil 

111 

555  3  55551 
;Ssus:t:i 

i::j5: 

itttn 

1  aaa* 
<  aaaa' 

1  Bill 

laaaai 

it::::::::;::::: 

:H::::::i 

aaaaa 

aaaaa 

aaaaaaaa* 

■IlliHH 

il 

aaaM  Saaa! 

lIHifiil 

i  mm 

ni 

ii 

5:5:5  »n! 

ill 

•■an 

•  aaai 

::::: 

•  ihi 

■aaaa* aaa aaaaa a* 

•  *****  aialiiH*! 

■  ■■■ 

■  ■■•a  aaaaa 
aaailMlat 

]>«•*>•••! 

csn; 

HE! 

::::::::: 

jim 

aaaaa  aaaaa 

It::!!::?1 

Uliiliu' 

::::::::: 

iiiii 

Ip  33 

aaiaSaaaai 

aaaaa  aaaa* 
aaaaa *a&ai 
aaaaa  *aaa« 

:u:ss 

laahaa 

aaaai 

iaapi 

•Bail 

aaaai 

*«*■! 

■  ••■•••••••Ukkl 

■aaaaa ■•**iia**i 
Uaaii  itaMMi** 
Haiti  l^kimiai 

3r:i 

:::::::::: 

«?:■ 

a- 

.ukkktki 

:  iiiii::! 

aaSai 

n:tf 

aTSailaaaaa 

:::»Ss»i 

liilili! 

•aaaaaaaa 

■■■*■■ *jj*a 
::::::::::: 

aaaa*  saaa* 
■aaa* kaaa* 
aaaaa  aaaa* 

:::::: 

>■•■44 

t:t ::: 

aaaai 

iiiii 

■  lakHi  laiikiiiMI 
>|..llll»l  **B  •• 

9. 

1  |  =  J 

aaaaa 

Hill 

••  *;;ss® 

:t*3iSiir 

*aaaa 

**■■• 

aaaa) 

::::::::: 
aaaaa  *aaa* 

::::::::: 

sisnssni 

mm 

nth::::: 

aaaa* aaaa 

iaaat aaa* 

l::l: 5:5: 

■saaa****** 

aaaa*  aaiaa 

3555 

::::: 

5555! 

ia*aaaaa*Bk*aaaii 
■■**•••■*»■ aaaaa 
:::::::::: 

M**aaia{i 
aaaaa  a'  «*• 

taaiaialaa 
aaaaa  aaaaa 

taSal 

IkiM 

::sk 

a«: 

::  ttjK: 

:: 

IHii 

:::::::::: 

:k::kk! 

nHlifiE 

n:sits:ra 

ttUiu::: 

5555; !55t! 

::::::::: 

Hijilni 

ktMHk***l 

:*:::: 

sis:::::::: 

kaaaa*  ■■Hi 

liSijtMul 

:::::::::: 

::::: 

SSiii 

k:» 

331 

■aaaaa ■••■•■■■•i 
UMiiiaaaiMkki 

i*>* a a •■**••■*•■ 
i ■**■■ aa**aaaa*. 

fiaaaa aaS{« 
.  .  aaaaa 

aaaaa aaaaa 

!ss::::i 

HtS 

m 

«::: 

aaaaa 

aaaaa 

aaaaa 

is:::::::! 

::::;k:;s 

aaaaa  aaaaa 

;sz:nstc 

aaaaa iaaaa 

tttnsnn 

aaaa*  aaaaa 
aaiaa •• *  » I 

Bkktkiutl 

:::: 

ill 

■■aaai ***** 

■  (MM  ■•**■ 

iliiillil 

aaaaa  aaaa ■ 

nsnt:::s 

■aaaa aaaaa 

::::::555: 

isSEes 

aa*ai 

aaaaa 

aaaaa 

iiiii 

aaiaa 

itaaai aa^ia aaaaa 

•■•II aBaata**** , 

'a*#** *■»*• •**••! 

:55:i:i::::ii:::J 

Sheie 

a aaaaaaaaa 

ihm:::: 

iHi 

;■•■■ 

aaaaa aaaaa 
aaaaa  aaaaa 

***;t**kki 

55::;  5553 

aaaia  aaai* 

;::!:55:5: 

::::::::::: 

533355: 

ttt:: 

aaaa* 

::::: 

laataaaaaaaaajlaa! 

la.E>»akll*ak|i 

Iill*:::::; 

nsnsn:; 

III:: 

*aaaa 

iasaa 

aaaaa 

aaaa*  aaaaa 
aaaaa aiaa* 
aaaaa aaaaa 

■aaaa  *a*ai 

liUilSjl! 

l;I:5:5:l:5 

SSSaaKS:* 

iiiii  iBiaa 

355! 

aaaai 

::::: 

•akii laaaiaSSSa 

■  ••••] 

PlSH-iH 

aaaaa  aaaSa 

::s:s:::» 

::::: 

aaaabaaaaa 

IkUilttaa 
l»aai  jaaai 

aaaaa 
aiaaa 
•  a»a> 

Caat*  aaaaa 

•aaaa aaaaa 
mmmmi aaaaa 

iill 

•ta*a*taaa*i 
■  aaaa* iasaa  i 

i :::::::::: 

ii 

ataaa 

::::: 

::::: 

aaaa* 

aaaa* 

***** 

***** 

*Saaip  »a  •**#•*:* 

aaaaitaua 
aaaaiiaan 
aaaaa aaaaa 
aaaaa  aaaaa 

:ss:: 

aaaaa  aaaaa 

saa>a aaaaa 

u:n 

aaaaa 

aaaaa 

yiiA 

n:::::n* 

■■■ii ■•••k 

aaaa* iaaaa 
Kaaa* aaaaa 
aaaaa  aaaaa 
aaaa*  aaaaa 

n:!i5:Hi 

liiiiiliiii! 

iaaat kaaaa 

aaaaa 
aaaaa 
■  •**■ 
■■•a* 
aaaa* 

•53! 

5:::i 

i 

SP^Si  Immt'i 
’**** 

■s-?  -v.^^K t  »■»■■ 
*•■»• 

B9Wi:eE; 

,***■ 

•* 

•-•»*>■* 

w*»S^Hp  •  •  »*»■ 

PfHgs&^^B*  aaaa* 

aaaaa  aaaaa 
•*»•• aaaaa 

aaaaa  aaaaa 

:::::::::: 

kaaaa aafiai 
aaaaa aaaaa 

aaaai »•*«» 

aaaaa lat*> 
iia*» aaaaa 

■HSHiSi 

:::t :::::: 

■•isiii-n 

3i:S 

::::: 

5::: 

::::: 

::2i 

l 

nth::::: 

nth:n:i 

aaaaa aaaaa 
aaaaa aaaaa 

Hippl 

jajja  Jaaai 

"tit 

ri 

::::: 

aaaaa 

aaaa* 

aaaaa 
::::: 
•  aaaa 

sin: 

nsnnsn 

nitssnn 

•ssnssni 

aaaaa  aaaa* 
Kaaaa  aaaaa 

-aa  aaaaa 
:::::::::: 

si: 

ataaa  aaaa* 

aaaa* aaaaa 
iaaaa  aaaaa 
aaaa* aaaaa 

aaia*  aaaia 

:n:::s:n 

nshstcn 

aaaaa aaaaa 

KKSMIK 

■aaaa  aaaaa 

nfi-Hg: 

:b::;sb 

aaa**  aaaa* 

EEiittE 

:::::::::: 

tnpnsn 

555:: 55555 
lit  55:553 

aaaaa  **a*i 
::::::::: 
Eaaa*  aSaai 
aaaa*  *a**i 

»•••' ‘ 

itiii 

IS 

•  ***a*  •••*■< 

|.s:C(tRt;! 

IlHHiliili 

iiUMliil 

! 11:55:53! 

i::::::::: 

;:::! 51555 

•aaa* aaaaa 
:::::::::: 

II 

mil 

•fkjj 

m 

il 

m 

aaaaa 

Iiiii 

CKH 

lOii 

> aaaaa  ■■•■■■•aa* 

t::EE::::«E 

•■■■••■a; »»*•*• 

11111 

ilppil 

***** 

IW»:t:: 

aaaaa  kaaaa 
kaaaa  kaaaa 

iitil 

i::H:::U 

aaaaa aaaa* 
aaaaa Saaa* 

ns|* 

aaaaa 

■  *lM*kk» 

snnnns 

aaaa*  aaaaa 

ntunn: 

553  55555 
53:55553 

5515 535! 

aaaia #**•< 
■  aaai  a*aai 

li*Mll**M| 

TOM 

iill 

thn 

USE 

s 

filPHttt::: 

#•*** 

■-*»• 

S^SVY. 

SiH' ::::: 
■Biiiiii 

::::i 

aaaaaaCaaa 

tsit: fill: 

aaaaa  aaaaa 
aaaaa  aaaaa 
aaaaa  akkak 

Hi 

aaaaa 

aaaaa 
«aa 
ka»  a 

iSaaa 

kaaaa 

kaaaa 

aaaaa 

aaaaa aaai i 

mWm 

ikii 

:::r. 

aaaa* 

aaaaa 

s 

iai 

1 

::: 

iS 

SiBS 

isntsi 

EiSH 

JJKiK 

::::::: 

aiaaaES 

tiMkli 

rat** 

ip  IB 

BKittil 

aaaa* iaaaa 
■Saaa aaiaa 
aaaa* ••*•* 

nnsnst: 

mllhti: 

iill 

mm 

li 

5535:55; 

aaaaa  aiaa: 

*****  a***a 

!fiaaS$BaaBi 

:h:t::::»i 

!:i:?rl‘illi 

Hi 

inn 

m 

n 

k;;; 

m 

335 

i::5: 

s 

Slily 
i  mm 

iillli 

s. 

:  i 

aaaaa  aaaaa 

iii::::5:; 

laaaft 

tstss 

•aaaaaaaa* 

iljUilST: 

!iiy 

Iasi 

::::: 

st::: 

txi?stht 
s;n:uv; 
hatitiu 
>»**■ •••«! 

tHEnnt 

KaKai iaaaa 
■••a* •••»• 
•  a  as* ••**■ 
••a* •*#** 
» • < • ■•«»« 

55:5:5535 

k:::H:k 

Bijiinit 

kk*ilU*»t 
kaaai *kk** 
**•«• ***** 

mm 

:nnn::j 

mm 

555!  55555  i 
:k;::u:ki 

ntaktsn: 

mm\ 

Ii 

pi 

oil! 

in 

BUS 

m 

:3::: 

KIR:  SHE 

»k*f • aaaaa 

iiiiii&i 

dll: 

laaSSaaZaS 

tisSmt! 

i:k; 

aaaaa 

ttn: 

kkkl •*•»• 
kaaaa aaaa* 

ninns:: 

'•***  t*». 

kaaaa  *aaaf 
»;m  ***** 

5135533 

KKt.::*: 

illi 

an 

ii 

11! 

I! 

Bl! 

«ir\  ;  fm 

iiiii Iiiii 

ns: 

mm 

m 

e::ei::s;: 

ESSL-HESt 

htnsnts 

::rr:::ss; 

::t:::nn 

irtEnH! 

t;uj 

H5!!ini! 

mm 

mm 

Pi 

P 

sps 

aXaaSaaaafi 
aaaaa  aaaaa 

•Azt:  it::: 
:;s|iK3 

ill 

itzi 

aaaa 
a»a  • 
aaaa 

tlkagUtaS 

mm 

Itiffli 

Uiii 

***** 

tUkl 

Ins: 

SSat  aaaa* 

hit::!::: 

iliiiii 

*a*k*  aaaa* 
■aaa*  aaaaa 
kaaa*  ***** 
kiaaa  aaaaa 

ntHstt: 

ssssssssts 

mm 

fEtijHg 

nt:n:3 

HI! 

iiom 

|l|j 

ip 

IS! 

ill 

Si 

W«*f 

n 

m 

■BiSi: 

m 

1 

mil 

***** 

•  ***• 

:xd 

Ip 

Mi 

SB 

mm 

mm 

raaffll 
555!  fuH! 

SI 

•TTT 
’  *  • 

iii 

■ 

lira  but; 
i:iS3!5!! 

laaSS  i 

sisti 

mm 

SHi 

St 

HlHli 

SHI 

iiiijj 

in 

Hill 

mm 

Si 

m 

0 

Flight  Con< 
Input  Show 
Ron-Integri 

lltio 
Lng  T 
sting 

n+TT 

n 

ha 

S 

S 

24  - 

t  th« 
ervo 

T7TTT 

Ai 

P 

ttttpttjt 

Pl« 

rornJ 

itoh 

.*J J --t 

10 

l,  D. 

’t  R« 
Rate 

rrrjTrrrj 

7 

epone 

Reei 

<K1 

1! 

)#  to 
>onse 

r  "  1 

FTTSTHTT 

r 

a  fit 
Droc 
0.0, 

T77TTTT 

ep  C 
pa  w 

Kv 

TTTT 

AC 

omt 

it) 

* 

TTTTf 

Tier 

1 

30, 

TT-.-prrrrrrr+T 

See, 

id 

the 

.4) 

87 


GA/EE/67-2 


|fe|{f#fe-fei:  :-i4lll  lijliii  lllli'  ij:;i|ilff|{  1|| 

iff  }  |.  Ji;i  ;iJ|.  : :  ‘  i  i .  i  l : !  i  . :  Pi-  !Ji;  !  ;•}  ft*!  ii  |[;  ft  i  :  ,  H  :?i  i  ; 

iff  tp|t|T ?:tr»Ti^  :-F  *r;:7?rrni  ffettitf  iljinril  dM  rriM 

i!i  illi  ilii  JIL  •i-!r.-:.-l;:niiil;i’ii:)i}  iliilil  iiil  Mil  ut;  III,  lit  lift  f  II  I  jit:  |J| 


■■  t ;  ;m  rp  jrTnrnrp?^  n 

B  mfll  m  g  ffS  I  -4*  if  { .  32ft 

Tnrinjrii^': -jlrjiru  ItsjiH  i:i!  ‘Hr  l^ilj^nTT  trlHtt 1- 
tUiiU  Hu  in  iliilil  iliini.jiiiiiikiLikil 


mm 


g  M 11  fepp?  HpipllMifpjJH'l' 

•!  tt::  ';,t  I’ri  ?■'**  :"t  'In  :!;)  ftti  :  f 

U4  -444  444!  flii  ;;4Uit4  444-4+44:.  414 iLi 44,04 lirrt  ifk  ulife  ft  Sii  Si  ill . 


t  rrn  mi 

M  ;!L 

ilj:  jiii  nil 

i til  Ii! : 

-44-  44*5: 

|  |tj|  |T;i  tijl 

I  -lil  j::i 

I  ll‘{  fej  i!ri 

•|i|  j-jj  ; 

j  i||l  jjH  jjjj 

;t.:  TH;  . 

♦  j  •  t  *  t  •  i  * 

1  till  lie  Iff? 

till 

jplkl  p 

i  til  i! ; !  1  • 

ii:i  1?[;  fit; 

pH'  V:  * 

f ji;  j!»|  •  • 

ji j  • : 

lif 

i!*i  i;|i  | 

'  tteinl  rill 

ill!  mi 

"  •  :  i*i  j 

j::;  ii*-  ! 

444-  |«  *»j»' ‘ I- 

77  ■  “!  i  ? 

:jij  ;  t;  ; 

rm  i  hi  j  M  1 


’  li 


to 

v\g,  i>.a 


ao  sec. 


Flight  Condition  y-\  —  Aircraft  Reanonse  to  a  Stop  Command 
Input  Showing*  That  the  pitch  Rate  Response  Remains  Steady 
with  the  Integrating  Servo  (Kf  »  10.0,  •  30,4) 
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Appendix  E 

Approximation  of  the  Aircraft  Pitch  Rate  Transfer  Funotion 
by  Finite  Difference  Equations 

Derivative  Approximations 

When  a  function  y(t)  is  tabulated  at  equal  intervals  (T  )  of 
an  independent  variable  t,  the  derivatives  of  the  function  at  any 
of  the  tabulation  points  may  be  obtained  by  the  use  of  finite  diff¬ 
erence  formulae.  These  formulas  all  involve  linear  combinations  of 
the  tabulated  values  in  the  neighborhood  of  the  point  where  a  deriv¬ 
ative  is  required  to  be  evaluated \  such  a  point  is  known  as  the 
pivotal  point.  In  general,  the  independent  variable  is  not  necess¬ 
arily  time,  nor  do  the  points  have  to  be  tabulated  at  equal  inter¬ 
vals*  However,  in  the  simulation  of  a  physical  system  by  means  of 
differential  equations  which  relate  the  time  variation  of  the  output 
to  the  input,  the  independent  variable  of  interest  is  time,  and  it 
is  a  praotioal  convenience  to  sample  the  system  behavior  at  equal 
increments  of  time. 

Differentiation  formulas  may  be  generated  by  starting  with  a 
Taylor  series  expansion  for  the  (usually  unknown)  function  in  the 
region  of  a  pivotal  point.  Then,  through  the  use  of  the  differen¬ 
tiation  operator  D  and  the  forward,  oentral,  backward,  and  averager 
operators  A ,  ft,  w,  and  respectively,  the  neoessary  weighting 
factors  for  the  tabulated  values  oan  be  found.  Details  of  the 
derivation  of  formulas  for  the  weighting  faotors  are  contained  in 
Ref.  19,  Ch.  ZI,  from  which  Table  IV  was  compiled. 

Construction  of  Difference  Equations 

As  an  example  in  the  use  of  Table  IV,  a  difference  equation 
will  be  constructed  for  the  airoraft  pitch  rate  transfer  function, 
equation  (C.4).  Re-writing  equation  (C.4)  in  differential  operator 
notation,  and  with  and  E ^  representing  and  respectively, 
at  the  pivotal  point  i, 

(D2  +  2^waD  +  w^)Yi  -  Mj(D>l/TR)Ei  (E.l) 
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Table  IV 

Weighting  of  Sample  Point  Values  for  Derivative 
Approximations  at  the  ith.  Point 


Order  of  Approx¬ 
imation  and  Type 

Scaled 

Derivative 

Sample  Points 

i-3  i-2  i-1  i  i+1  i+2 

1st.  Forward 

T*t> 

-1  1 

T.V 

1-2  1 

1st.  Central 

ar4D 

-10  1 

r/D* 

1-2  1 

1st,  Backward 

r,D 

1 

M 

r.V 

1-2  1 

2nd.  Backward 

2rst> 

1  -4  3 

r/D* 

-1  4-52 

To  approximate  equation  (E,l)  by  means  of  central  differences, 
it  is  first  multiplied  on  the  left  hand  side  (LHS)  and  the  right 
hand  side  (RH3)  by  Tg»  Then,  from  Table  IV, 

D^Y  ■  (Y  -2Y+Y  ) 

's  i  '  *i-l  i  i+1' 


"Var6  ^~Yi-l  +  Yi+1^ 


and, 


jy,  ft. 

T_  w 


2Y< 
a  i 


2^,2  y 

Vs  *i 


Therefore,  by  adding  the  above  equations,  and  oolleoting  like 
terms, 

ue  .  a  -  v.Wx  *  <-X-  2>ri  *  <i 


Similarly, 

HHS  •  (Ei+1  -  E^)  +  ^  *t 


Thus,  equating  LHS  and  RKS,  and  solving  for  Y^ 
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♦ 


let#  Order  Central  Differenoee 

h.l  ■  [  (?  -  >  T1  -  (1  -  >  Vl 

>i'i,  rJn  ] 

♦  ~f~  (Em  ‘  Ei-i>  ♦  ~ t  Ei|  <E-2) 

a  J 

In  a  similar,  way,  by  using  alternative  approximations  for  the 
derivatives  from  Table  XV*  the  following  equations  may  be  obtained* 

. '  i  i-  '  . 

1st*  Order  Forward  Differences 


m  ’  (24*.'r.  -  -  1  )Yi-i+  )t. 


2 

♦  r,“(<Ei  -  Vi>  *  Ei-i  <E-*> 

*  a 


1st.  Order  Backward  Differences 
T,  - 


-  . 'V - ~TT  [  <2  +  2*aWaV>Yi“  Vi 

1  +  2  *a1,»‘rB  +  ars  *>  a  a  s  i  il 

rf2  U  1 

♦  WW  -V  *  -If  E1+1J  <*.4> 

2nd.  Order  Backward  Differences 

Yifl  *  - - - - TTJS  [  (10+8^w  T  )Y  +  2Y 

4  +  6^war#  +  2wa1f8  l  ■»  a  s  i  i-2 

-(8+2*u\r9)r lml  -  4*^ 


2rth  1 

•“T-->«i+ij  Os. 


5) 


An  additional  equation  may  be  obtained  by  using  a  let.  order 
backward  difference  approximation  for  the  eeoond  derivative,  in 
conjunction  with  2nd.  order  backward  differences  for  the  first 
derivative.  This  mixed  order  backward  difference  approximation 
warn  used  by  Stein  (Ref.  21i66). 


| 


i 
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Mixed  Order  Backward  Differences 


ifl 


,  77— f(4t8i.%n)Yi- 
2+6v.n+  2,aT.  1 


4“lr.Ei 


-  (2  +  2y  w  Y  )  Y4  -1  +  M.  T  E.  _ 
'  '’a  a  *s  i-1  *  s  i-1 


(3  Vl 


(E.6) 


Iteration  of  Difference  Equations 

The  forcing  function  in  equations  (E.2)  through  (E.6)  is 
represented  by  the  sampled  values  (n  *  ...  -1,  0,  1,  ...)  at 

integral  multiples  of  T  .  At  these  sampling  times,  as  each  new 

B 

value  of  En  becomes  available,  an  iteration  may  be  made  to  find  the 

computed  value  of  Y^.  Thus,  when  the  subscript  i  is  set  to  zero 

in  the  above  equations,  expressions  for  the  computed  values  of 

are  obtained.  If  the  equations  are  iterated  for  zero  initial 

conditions  all  the  E  and  Y  for  n<0  are  zero,  so  that  Y,  can 

immediately  be  calculated.  However,  for  non-zero  initial  conditions 

two  or  three  values  (depending  on  which  equation  is  ussd)  of  end 

Y  for  n<  0  are  required  to  be  known  before  the  iteration  can  be 
n 

started.  In  the  closed  loop  simulation  program  for  the  complete 
oontrol  system  (Appendix  A)  it  was  therefore  necessary  to  use 
temporarily  stored  values  of  elevator  deflection  angle  and  aircraft 
pitch  rate  (from  the  rate  gyro)  in  order  to  start  the  difference 
equations. 

Aoouraoy  of  Difference  Equations 

To  test  whioh  of  the  difference  equation  approximations  to  the 
airoraft  transfer  function  was  the  most  aoourate  a  oomputer  program 
was  written  for  the  AFIT  IBM  1620  maohine.  In  this  program  the 
airoraft  pitoh  rate,  oaloulsted  by  one  of  the  differenoe  equations 
(E.2)  -  (E.6),  was  oompared  with  an  analytioally  oaloulated  airoraft 
response.  A  listing  of  the  program  statements  is  given  at  the  end 
of  this  Appendix,  and  an  outline  of  the  oomputatione  made  Is  shown 
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in  Fig.  E.l. 

Since  this  testing  was  done  with  no  dosed  loop  control  system 
around  the  aircraft,  it  was  necessary  to  assume  an  arbitrary  elev¬ 
ator  motion.  Equation  (E.7)  in  Table  V  was  ohosen  to  define  the 
elevator  deflection  as  the  response  of  a  2nd.  order  transfer 
function  (damping  factor  J'g,  natural  frequency  w@)  to  a  step  input 
of  amplitude  S.  An  analytical  function  for  the  aircraft  response, 
equation  (E.8),  to  this  elevator  motion  was  obtained  with  the  aid 
of  Laplace  transform  tables  (Ref.  13*87).  Plots  of  equations  (E.7) 
and  (E.8)  are  contained  in  Fig.  E.2,  together  with  the  values  of 
parameters  which  were  used  to  generate  them.  It  may  be  noted  that 
the  values  chosen  for  the  aircraft  parameters  were  representative 
of  a  typical  Flight  Condition,  so  that  they  lie  towards  the  centers 
of  the  ranges  of  variation  shown  in  Figs.  2,3,  and  4*  Also,  by 
comparing  the  elevator  motion  for  the  first  0.5  sec.  in  Fig.  E.2 
with  that  in  Fig.  D.8 ,  it  is  seen  that  the  input  function  used  to 
test  the  difference  equations  was  a  realistio  approximation  to 
elevator  motion  whioh  could  occur  in  the  closed  loop  system.  Hence, 
the  results  of  testing  the  accuracy  of  the  difference  equations 
under  open  loop  conditions  were  expected  to  be  equally  valid  for 
the  closed  loop. 

Figures  E.3  and  E. 4  show  the  normalized  r.m.s.  errors  of  the 
difference  equations  as  a  function  of  time.  The  apparently  large 
initial  errors  are  due  to  the  faot  that  the  normalization  was  made 
with  respect  to  the  aircraft  pitoh  rate  (as  indicated  in  Fig.  E.l), 
whioh  was  zero  at  time  zero.  From  these  two  Figures  it  is  seen  that 
1st.  order  central  differences  give  more  accurate  results  than  any 
of  the  other  1st.  order  approximations,  and  are  even  more  aoourate 
than  2nd.  order  baokward  differences.  The  curve  for  1st.  order 
backward  differences  followed  that  for  1st,,  order  forward  differences 
in  Fig.  E.3,  and  was  therefore  not  plotted. 

A  particular  feature  in  Figs.  E.3  and  E.4  is  that  the  errors 
for  all  but  the  central  differences  start  to  inorease  when  the  first 
derivative  of  the  airoraft  pitoh  rate  (Fig.  E.2)  ohanges  sign.  This 
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is  an  indication  that  the  approximation  for  the  seoond  derivative 

by  central  differences  (Table  IV)  is  the  most  aoourate,  a  result 

wnich  may  be  derived  theoretically  (Ref.  19»86).  For  x  “  0.01  see. 

8 

the  normalised  r.m.a.  error  of  the  central  difference  equation  was 
found  to  be  0.01  per  cent  after  ?.0  sec.  (i.e.  after  200  iterations); 
the  actual  error  at  this  time  was  0*0008  deg./seo.,  when  the 
anal/tica!  value  of  pitch  rate  (^ggo)  was  ”1*9850  deg./sec.,  so 
that  the  terminal  error  was  0.04  per  cent. 
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Table  V 

Analytical  Expressions  for  Elevator  Deflection  and  Aircraft 
Pitch  Rate  Used  for  Testing  the  Accuracy  of 
Difference  Equation  Solutions 


Elevator  Deflection 


i9(t)  s  +  (B«7) 

u»Aer«, 

A  .  — L-. 

8,  «  w, 

Cj  *  “W/-.S* 

o,  »  f.n- ' 

Aircraft  Pitch  Rate 


tjt)  «  - Siafc.t-O,)  +^C  atAii»(C*e-dO]  (E.8) 

w>A«t«, 

A  .  _ «.u;r.1]1'* _ 

[«-***«•*$,  *4,  *a  ]  ** 

fl,  *•!*«•»• 

C,  •  W* 

A  «  tan'  JhH  -  4an‘  .1  ay  (3a0)« -  3« _ 

+  *Z}al  -u£ 

-w*  *\*g 

The  oonstante  Ag,  Bg,  Cg,  and  Dg  are  obtained  by  interchanging 
the  subscripts  a  and  s  in  the  expressions  for  A^»  B^f  C^»  and  D^» 
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pRts&RAM  INPUT  OATA 

i,  fteflttseMraT/ve  v*i.uts  M  tm*  mm cmmmr 
mm hah e nits  i  Ti,  5*  »  **>- 


.  c**«r-*wr*  r*  %y#rHKtix.€  A  A  o4#«4 
ZLLiRToK  *£*P*M**\  &, 

1 3.  srmrlk  4w*  trtMTft/  iMT*AY*L  :  Tf 


AiRCRAFT 

6«(*)  •* 
e*uar/*w  (e.9) 


OIFFEREMCE 

EQUATION 

Ml  M*  Kt«U- 

(c.a)  -  (€•*) 


44*41* 


J! - - 

ELEVATOR 

S,(fc)  #i*«*  *y 

Csusrienf  (1.7) 

4*  4 
rawer***  «r 
rmt 


*•4*4  i./«*S 
TgJ.  £44#£ 

£(Vfl s>* 


!{■ 


Fig.  K.l 


Calculations  Mads  in  ths  Program  Ussd  to  Tsst  ths  Aoouraoy 
of  Differsnos  Equation  Approximations  to  ths 
Aircraft  Transfer  Punotion 
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((Uj./sec.) 


Fig.  E.2 

Elevator  Deflection  and  Analytically  Calculated  Aircraft 
Response  for  Repreeentative  Values  of  Aircraft  Parameters 
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A/<rrma(«x«dL 

r.m.s. 

%error 


/•o  hr  2-0  Me. 

Iteration  Length 


Fig.  E.3 


Reiatiye  Aocuraoy  of  Various  Difference  Equation  Approximations 
Showing  that  Central  Differences  are  the  Most  Acourate 
Even  #hen  the  Iteration  Length  is  Longest 
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Msreiet<*eJl  O' S' 

r.ro.s. 

X error 


ir  ao  see. 

Iteration  Length 


Fig*  E.  4 


Aoouracy  of  let*  Order  Central  Differences  Compared  fith 
2nd.  Order  Backward  Differences  Showing  that  the 
Central  Differences  are  More  Aoourate 
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C  ACCURACY  OF  DIFFERENCE  EQUATION  APPROXIMATIONS 


DIMENSION  S I GH ( 2 ) 

READ  *  ELDEG*  WEL.  ZETEL 
READ ♦  H 

C  REPRESENTATIVE  AIRCRAFT  CONSTANTS 

DELTM«10.0 
WAC»3.0 
ZETAC*0.1414 
TAC*5.0 

C  CONSTANTS  FOR  THE  ACTUAL  PITCH  EQUATION.  Al-D2.<8) 

S* ( ELDEG  *DELTM ) / <  TAC«WAC*WAC ) 

W1*WAC 

W2*WEL 

Z1*ZETAC 

Z2-ZETEL 

CALL  SIZE 

GO  TO  13 

SUBROUTINE  SIZE 

F1«W1*Z1 

F2»W2*Z2 

HNUM«W2*W2*SQRT( 1.0-2.0*TAC*Fl+TAC*TAC*Wl*Wl > 
RT«SQRT( 1.0-Z1*Z1) 

C*W1*RT 

D»W2*W2-2.0*Fl*F2+2.0*Fl*Fl-Wl*Wl 
DEN»PT*SQRT<D*D+4.0*C*C*IFl-F2)*tFl-F2) ) 

A  •  HNUM/DEN 
RETURN 

13  Al-A 
Bl-Fl 
C1»C 
I«1 

11  TOP-RT 

BOT*-Zl 
CALL  ANGLE 
GO  TO  14 

SUBROUTINE  ANGLE 
IF(B0T)4  *3*5 

4  IF(TOP)9*3*7 

9  P*3. 1415926+ATAN ( ABS ( TOP/BOT ) ) 

GO  TO  2 

7  P*3. 1415926-ATAN ( ABS  I TOP/BOT I ) 

GO  TO  2 

3  IFITOP ) 8  *3  *6 

8  P*6. 2831853-ATAN ( ABS (TOP/BOT ) ) 

GO  TO  2 

6  P-ATAN  <  ABS I TOP/BOT  J ) 

GO  TO  2 

3  PUNCH  50 

50  FORMAT (UHANGLE  ERROR) 

2  RETURN 
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14  P1«P 

T0P*TAC*C 
BOT*1*0-TAC*F1 
CALL  ANGLE 
P2-P 

TOP“2*0*C* (F1-F2 ) 

BOT-D 

CALL  ANGLE 

SIGHm»Pl-P2-P 

I*I+1 

IF( 1-3) 10* 12  *  12 
10  Wl-WEL 
W2-WAC 
Zl-ZETEL 
Z2*ZETAC 
CALL  SIZE 
A2»A 
B2*F1 
C2-C 
GO  TO  11 
12  D1*SIGH( 1 ) 

D2«SIGH( 2 ) 

CONSTANTS  FOR  ELEVATOR  DEFLECTION  EQUATION*  A3~D3*(4) 
A3«1*0/RT 
B3-F1 
C3*C 
D3*P1 

SHORT-PERIOD  CONSTANTS 
ZWAY*ZETAC*WAC 
RTAY-l.O/TAC 
WAY* W AC 

COEFFICIENTS  FOR  CENTRAL  DIFFERENCE  EQUATION 
CY*1.0/(1*0+H»ZWAY) 

CYE31*H#DMY#0*  3#CY 
CYE2«H*HPDNY»CY#RTAY 
CY2*C Y* ( 2  *  0-H*H#WAY»WAV  > 

CY1*( 1*0-H*ZWAY)*CY 

INITIAL  CONDITIONS 
PUNCH  96 

36  F0RNATO3H  TINE  ELEV*  AC*P*R*  CONP*P*R*« 

C  18H  ERROR  RMSY  P.C.*//) 

AC3«0*0 

AC2-0.0 

£3*0*0 

£2*0*0 

Y3-0.0 

Y2*0*0 

SNSQA*0*0 

SNSQY*0*0 
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ITERATION  OF  THE  DIFFERENCE  EQUATION 
T*0.  0 

L“0  .  'V;  ■  '!  . 

61  L-L+l  '  ;■  Vi'. 

90  T-T+H  V 

AC1*AC2 
AC2* AC3 

APT1«A1*EXP<-B1*T>*SINK1*T*D1) 

APT2*A2*EXP<-B2*T>*SIN«C2*T-02» 

AC3«S*f  1.0+APT1+APT2) 

E1-E2  ■ 

E2-E3 

£3*ELDE<3«(UQ+A3*EXP<-83*Tr*£I>4<C3*T-D3n 

Y1-Y2 

Y2-Y3 

Y3*CYE31*<E3-E1  i-KYE2*E2*»-CY2*Y2-CYl*Yl 
SMSQA°$MSQA+AC3*AC3 
SMSQY»SM$QY+<Y3-AC3 >#IY3-AC3I 
CNT«l 

TEST  «ABS<  T-0. 1*CNT> 

IF ( TEST-0.00001 >32»32#90 
32  RMSYP*100.0*S0RTCSMS0Y/SMSQA) 

ERR»Y3-AC3 

PUNCH  34.  T.  E3»  AC3.  Y3»  ERR »  RMSYP 
34  FORMAT  t  F3.2*F7«2  »F9*4. F20.4.2F9.4) 

IFCL-20 161.99.99 
99  STOP 
END 
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